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Heat transfer to air flowing through packed tubes 
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Chemical Engineering Laboratory, Applied Chemistry Department, College of Technology, Manchester, England 

Abstract—Various methods proposed for the expression of heat transfer rates between the 
containing wall of a packed tube and a fluid flowing through the packing have been reviewed. 
By using a simplified method, additional experimental data have been obtained which indicate 
that the characteristics of a packed bed may be adequately described in terms of a uniform 
effective conductivity and an effective wall coefficient, both of which are essentially proportional 
to the fluid velocity. The variation of the local overall heat transfer coefficeint with positions 
along the length of the bed has been examined. 


Résumé 
déchange thermique entre les parois dune colonne a garnissage et le fluide qui sy écoule. L’emploi 


Les auteurs ont repris les différentes méthodes proposées pour lobtention des vitesses 


dune méthode simplifiée leur a permis dobtenir des données expérimentales supplémentaires 
prouvant que les caracteristiqués d’un lit a garnissage sont fonction d'une conductibilité effective 
uniforme et d'un coefficient de paroi effectif : cex deux termes sont, en fait, proportionnels a la 
vitesse du fluide. Les auteurs ont examiné la variation locale du coefficient d’échange thermique 


global en fonction des positions le long du lit. 


Later workers [11, 15, 
uniform effective conductivity for the bed derived 


Errorts to provide a correlation for the rate of 21] preferred to use a 
heat transfer from the tube wall to a fluid flowing 
through a packed tube have shown a progressive from overall temperature differences with the 
Karly workers [5, 7, 13, 14] ex- 


pressed their results as a uniform effective overall 


assumption that the velocity of the percolating 
The 


two methods of expressing heat transfer data can 


developemnt. 
fluid was uniform over the bed cross-section. 
heat transfer coetlicient h, analogous to that for 
be related, but are both limited in usefulness by 
Smiru et al, [4, 9, 18] derived 
values of the local effective conductivities of the 


a fluid moving through an empty tube, namely 


Q 


their definitions. 


packed volume by differentiation of experimental 
the 
Their 


suitable mean temperature distributions, which 


the tube wall and the fluid 


is a temperature express 


AO, 


between 


where 
difference conditions in the interior of the packing. 
stream. A full representation has not been reported data should be comparable with effective diffusi- 


for the effect on this coeflicient of variations of — vities in packed beds obtained by other methods, 


the several parameters used to describe the bed 
and the fluid. Whilst there is no essential in- 
accuracy introduced by the use of a coeflicient 
defined in this way, it would probably be too 
complex a function of operating variables for 


normal use. 


* Present address : 


An effective conductivity of the packed volume 
alone, Le. excluding effects due to the presence 
of the wall, should give a better expression of the 
properties of the packing and be independent of 
the geometry of the wall surface. A logical step 


appears to be the separation of the effects of the 
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boundary wall to provide an expression of the 
transfer properties of the packed volume in 
terms of the characteristics of its components, 
but it must be realized that the wall’s effects ts 
not limited to its effect on the heat transfer 
characteristics in its near vicinity. The fluid 
velocity distribution throughout the bed ts 
modified by the conditions near the wall [19], se 
that the effective conductivity of the packed 
volume, k,, being dependent on this velocity 
distribution, Is also changed. The simple sub- 
division of the system into a local wall coetlicient, 
h,,. and a packing characteristic, &,, must there- 
fore be considered ais only a second approximation 
to the true conditions. The earliest attempts to 
make this subdivision by Coperty and Marsa. 
and by Brorz | 3) can be criticized in matters 
of detail in that entrance effects were neglected, 
uniform fluid velocity was assumed, and simplify- 
ing assumptions were made about the form of the 
temperature distribution ; but the method 
represents a definite advance. The present 
analyses employ essentially the same bases of 
expression, but use a simpler analy tical procedure. 

Earlier work [5] appeared to prove that the 
overall heat transfer coeflicient, A, \ iried with 
the depth of the bed. Although this had not 
heen contirmed by other reports it seemed 
reasonable to expect that some change in / must 
take place in the inlet sections of the bed while 
changes are occurring mm the velocity and tem 
perature profiles in the fluid. In contrast &, for 
the packed space should show smaller variations 
resulting from changes in the velocity profil 
only as the fluid enters the bed. Any subseque nt 
changes in velocity distribution due to inhc mo 
geneities of the packing or at the exit from the 
hed should also be reflected in further changes 
in Ad and toa smaller extent in k.. These variations 
within a given bed cannot be assessed correctly 
by a study of the behaviour of a series ol beds of 
differing length, but should be obtained from 
measurements of local conditions within a singk 
hed. The present work included an attempt to 


make these measurements. 


EXPERIMENTAL Bases 


The experimental investigation of both the 
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variation of h with position along the length of 
the packed bed and the subdivision of the heat 
transfer characteristics of the system into /,, and 
k, are greatly simplified by imposing the con- 
dition of constant heat flux per unit area of wall 
surface rather than that of constant temperature 
at the wall. The following advantages are 
obtained : 

(i) The variation of surface temperature along 
the length of the tube wall is easier to measure than 
either the rise in temperature of the fluid stream 


or the local heat flux from the wall. 


(11) Longer lengths of packed bed may con- 
vermently he studied and reasonable radial and 
longitudinal temperature distributions  main- 


tained. 


(iii) The radial temperature distribution can be 
shown to be parabolic if the uniformity of fluid 
velocity is accepted, and this fact offers a con 
venient method of graphical expression and 
averaging. The need for graphical differentiation 

of experimental curves to obtain is 

or or 
removed with consequent gains in speed and 
accuracy. The assumption by Brorz [3] that the 
radial te mperature distribution ts parabolic when 
heating is maintained with constant wall tem- 
perature ts ineorrect. The extent of the error in 
this assumption varies with the value of the 
Graetz number, cW kL, for the thud. 


rradient, 


(iv) The longitudinal te mperature 
dd 


is a constant, independent of radius when & ts 
large. This quantity can therefore be obtained 
more accurate ly than for the case of heat transfer 
with constant wall temperature. In the latter case 


(0, 0)| 


the value of Is required, which 


is theoretically most accurately obtaimed when 2 
is large, but (@, is. then increasingly 


difficult to measure with the requisite precision, 


(v) Hitherto, although the theoretical tempera- 


ture distribution was available for the case of 


constant wall temperature, the case for constant 
heat flux at the wall had not been analysed, The 
late Dr. S. R. Tress [20] of this College (¢ ollege 
of Technology, Manchester) provided the solution 
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for the case of fluid with uniform thermal con- 
ductivity with uniform temperature at entrance, 
and constant velocity v at all positions of r and 2, 
the tube wall at radius a being heated to give a 
constant heat flux qg’ per unit area of tube surface 
at all wall positions. The temperature @ at any 


position within the tube at radius r and at 


(a) 


Fic. 1. 


distance x from the inlet to the heated section 
is given by 


2gor ,qr” qa 


= 4k 


(A, \? 


Bessel function of zero order of first 
kind, 
Bessel function of first 
kind, 
A nth root of Bessel function J (A) 0. 


n 


avk 


J, (A,) 


order of first 


EXPERIMENTAL METHOD 

The apparatus in Fig. 1 consisted essentially of a 
tube A made from nickel foil of 534 in. length and 
1% in. diameter with 0-002 in. wall thickness. The 
tube was packed with glass spheres Q, of mean 
The 


was 


diameter 0-173 in. and bed porosity 0-37. 


controlled mid-section of 29}in. length 
heated by the passage of a direct electric current 
through the tube material comprising the section. 
The tube was effectively lagged and strengthened 
by an annular jacket of 0-45 in. paraflin wax B 
and an outer jacket of 1 in. kieselguhr powder C. 
This lagging was enclosed in an air jacket, formed 


Finally 


from cardboard tube N with spiral P. 


Sf 


Heated tube apparatus. 


247 


n=@ 
2 
= 
| 
OL, 
5 
956 
c 8 
® 
| 
‘ A 8 Cc 
a 
ASO 
\ 4 \ / e 
| | | 
= 


J. H. Qurixton and J. ANDERSON STORROW 


jacket 
winding D 


this 
electric 


through was 


the air temperature 


carefully controlled by an 


placed on the outside surface of the kieselguhr 


lagging. the rate of heat input and the air ve locity 
being adjusted to make the jacket temperature as 
nearly equal as possible to that of the tube wall 
at all positions along the tube length. Heat losses 
were thereby reduced to a minimum . 

Tube 
three 
thin-walled glass quills F strapped to the outside 
of the foil tube, 


circumference. 


wall temperatures were measured by 


movable thermocouples k cm losed in 
and equally spaced irround the 


The 


selected to be of uniform bore 


quills were cure fully 


just large enough 
the the 


thermoce uple were 


movement of rmocouple 
beads within them. The 
made from Nichrome (80 Ni 20 Cr) and 
40 Ni) 


fully fused together to a sound globular junction 


to allow the 


Advance 


(oo Cu wires of 0-005 In. diameter care- 
of about 0-02 in. diameter to give a reat fit in 
the appropriate quill. From a study of the relative 
and its 
that the 


tight-fitting beads held by wires in tension would 


heat exchange rates between the bead 


immediate surroundings, it was felt 


allow an accurate mexsurement of local tem- 


perature, The fine quills strapped to the foul 
tubs 


take up the temperature of the insula 


insulating wax would 


tion at that 


and embedded in the 
position. At a point so near to the tube wall, 
when the radial temperature gradient was rrnaacle 
small bw adjusting the jacket to verv small heat 
loss into the lagging, the thermocouple tempera 
ture could be taken as that of the tube wall. With 
the quill arrangement and moving thermocouples 
a continuous measurment of tube w 


tub 


tempera- 
tures was obtained along the length, giving 


detailed and 
obtained 


a far accurate distribution 
than 


thermocouples at 


nore 


could be with many separate 


fixed positions around and 
along the tube length. Similar thermocouples E,, 
face of 
the pari flin wax to measure the temperature O, 
at the kiesel- 


guhr. 


were inserted in quills F’ on the outs 


interface between the wax and the 

The measurement of outlet air temperature was 
found to be critical. <A 
upward flow of air was the indirect measurement 
of the 


very method used for 


mean temperature from experimentally 
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determined profiles using a thermocouple which 
could traverse the top face of the bed in all 
directions. This probe thermocouple is shown in 
Fig. 1 and consisted of a brass tube J supported 
by a universal joint K, carrying at one end the 
thermocouple junction M supported at the tip 
of a thin-drawn-out twin-bore silica quill L. The 
position of M was measured on a circular scale 8 
above the porter at the top end of J. 

As the heat loss could not be reduced precisely 
to zero at all points on the tube wall at one time, 
a small correction was applied to the heat input 
at each position of measurement along the tube. 
The apparatus was calibrated for heat loss by 
experiments designed to provide a known heat 
through the 
temperature distribution similar to the conditions 
The 


temperature difference between the tube wall at 


flux insulation and with an axial 


in which the apparatus would operate. 
Oy and the lagging at @,; was noted for each 
heat flux per unit area at each position along the 
tube length, and it was found possible to express 
the heat loss for each position as a simple linear 
function of the t mperature difference. Exce pt for 
positions within an inch of the ends of the heated 
section, where mean values had to be used, the 


results were consistent as shown in Fig. 2. The 


;. 2. Heat loss calibration chart. 
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data for a selected position fall on a straight line 
within 0-02 B.Th.U. /in.hr. The uncertainty 
introduced at the ends is exemplified by the data 
for a hin. and « = 29}in. The abscissa in 
Fig. 2 is the heat flux per unit time per unit 
length of tube at the stated axial position x along 
the tube. 


experimental range of 10 to 50°C a constant 


For the thermocouples used in the 


conversion factor of 22-4°C ‘millivolt could be 
used with a maximum error of 0.4% in derived 
heat transfer quantities. The total heat loss was 
always reduced to less than 5% of the total heat 
input from the foil, which varied from 10 to 
130 B.Th.U. /hr.. and the estimate of the total 
loss was considered correct to within 0-5 to 1-0 
B.Th.U. br. over the same range. The heat flux 
at each position along the tube was obtained by 
subtracting the rate of heat loss from the rate 
of heat generation as derived from the wattage 
dissipated in the foil. The corresponding tem- 
perature of the air at each position was calculated 
from a heat balance equating the total heat input 
from the inlet to position a to the heat absorbed 
by the air stream 
G (Om: — 70" = 
0 

where gq, is the heat loss per unit length of the 
tube and varies with position r. When the overall 
heat balance failed to meet the measured inlet 
and outlet temperatures were always assumed 
correct and the heat capacity of the air (G-c) 
adjusted accordingly. A few tests in which the 
heat balance discrepancy reached 10% were 


discarded. 


EXPERIMENTAL RESULTS 


Before beginning the present investigations on 
a packed tube, the apparatus was tested without 
the packing, and coefficients were derived for 
fluid moving through an empty tube. The results 
were found to agree with existing experimental 
data and theorteical predictions. The expert- 
mental results for a ty pical tests on a packed tube 
are shown in Fig. 3. The data in the top section 
show the thermocouple M readings for traverses 
across four equally spaced radial positions. The 
lower section of Fig. 3 shows the distribution of 


mv 
z 
= 1-24 
83 
1-32 
2% MEAN OUTLET = g 
zi 
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z 
as 
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of 5 10 270 30 


OISTANCE FROM INLET 


Fic. 3. Measured temperature distributions across and 
along the heated tube. 


wall and lagging temperatures at positions along 
the tube, each te mperature being measured by the 
three thermocouples equally spaced around the 
same circumference. Within the heated section 
shown by the vertical lines the agreement for 
temperatures Is good and the controlled lagging 
has made the lagging temperature slightly higher 
than the wall temperature. As shown in the 
exemplifying calculation in Table 1, the heat 
input was small, but the sensitivity of the 
apparatus is clear from the rapid changes of 
temperature distribution at the ends of the 
heating section of the tube. The mean tempera- 
ture of the air stream is shown by the ringed 
points for the measured values at the inlet and 
outlet to the bed. 

The variation in overall coefficient A with axial 
distance from the inlet to the heated section Is 
shown on Fig. 4. The experimental results cover 
a wide range of flow rates. In all cases the only 


significant variations occurred in the first ten 
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Table 1. Run No. 256. Section A : Measurements along tube length. 


Distances from inlet Temperatures Effective heat 
— —— —_—— Heat Heat Total Air transfer 
Section Section gain absorbed heat Tempera- AO coefficient 
mid- outlet Ox O, absorbed ture (h) 
position position B. B.Th.l ‘) B. B.TAU. ) 
(in.) (in.) (mV) (mV) (mV) (nV) ft? hr. °F. 


hr. 


29} 294 
284 


67-38 195 144 12-3 
66-13 179 119 143 
63-69 ‘157 115 14-7 
241 61-29 135 106 15-5 


1-374 
1- 
1- 
1 
223 1-245 58-04 115 108 14:3 
1 
1- 


358 
335 


339 
298 


sit 
— 
— 


26} 206 
254 
234 1-177 
204 1-126 


l- 

274 1-272 

1- 

1 


179 46 56-73 O76 102 148 
130 50-27 109 14-0 
174 1-066 O97 re 43:74 0-959 108 142 
14) 1011 027 87-15 0-900 112 13:8 
11} 0-962 0-990 P86 30-53 0-840 123 13-0 
8} 0-898 0-926 23-67 0-776 122 13-4 

54 7 0-843 0-876 16-76 0-713 130 12-9 

23 0-774 O-S14 9-50 0-648 126 13-2 

0-725 0-775 2-37 0-606 119 13-9 

0 0 0-595 


+) 


Note—Nichrome-Advance thermocouples 40-4°F /mV. 


Column 1 : Distance from air inlet of heated section to measuring position for temperatures in columns 3, 4. 

Column 2 : The outlet position of each section for which column 1 represents the mid-position. 

Columns 3,4: Mean values in each case from each set of three thermocouples for lagging and wall temperature. 

Column 5 : Heat gain over the tube section between column 2 position and the position below it; each being 
assessed from the loss or gain rate from the graph (Fig. 2) appropriate to column 1 position, 

Column 6 : Corrected heat absorbed for the tube section in column 5 based on a constant heat input of 2-177 B.Th.U. 
per hr. per inch of tube length deduced from input electrical power. 

Column 7 : Total corrected heat input from inlet to position in column 2. 


Column 8 : Air temperature at column 1 position derived from column 7 followed by interpolation to column 1 
position. The temperature is given as the equivalent millivoltage for the experimental thermocouples 
with cold junctions in pure ice. 

Column 9 : The effective local temperature difference A@ (a 
microvolts for the thermocouples. 


« — OG») at mid-positions in column 1, expressed as 


Column 10; ~— The local transfer coefficient A from the heat transferred to the air round the mid-position in column 1, 


Section B: Radial measurements with moving outlet thermocouple. 


Temperatures are given as thermocouple millivoltage with the cold junction in pure ice. The temperatures were 
measured at eight equally spaced radial positions across the radius of 13/16 in. Each radial position n here and in 


13 \. 
Fig. 3 (top section) thus corresponds to a position {n Ss jin from the tube wall. The measurements were taken 


over four radii equally spaced round the tube section ; A, B and C, D being pairs on the same diameter. 


Radial position 0 3 5 6 7 

1-154 1-205 1-210 
1-180 1-232 1-250 
1-165 ‘ 1-206 1-210 
1-164 1/183 


Temperatures 


q 
— 
| A 
| 
VOL. 
1956 
q 
fat 
q 
c 
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1-230 
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Section C: Derived results 
Total heat input from foil = 64-2 B.Th.U. /hr. 
Heat gain by air from lagging = 3-17 B.Th.U. /hr. 
Inlet air temperature = 0-595 mV (57-2°F). 
Outlet air temperature 1-200 mV. (81-6°F). 


Heat absorbed from air flow rate and temperature rise = 63-2 B.Th.U, /hr. 
0 
5, = 865°F /ft.? ; 9-8°F /ft. 
or 
‘= Gd 
2340 ; | 250 
7 G = 738 Ib. /ft.*hr. 
k, = 0-50 B.Th.U.ft. /ft.*hr.°F. 
h,, = 26-5 B.Th.U, /ft.*hr.°F, 
AD 120 
k 


inches, i.e. six tube diameters, and depended on temperature drop (@,,— @,) leading to discon- 
flow rate. Only minor variations existed in the  tinuity in the temperature distribution in the 
remaining length of the tube. idealized model. From a heat balance on the 
annulus of thickness dr at radius r and length 8x 

ANALYSIS FOR Errective Conpvuctivitry 
AND Wau COEFFICIENT 


. 
br da (2) 


* 
or oa 


or 


The temperature distribution shown in Fig. 5 is 
consistent with the conception of a uniform After the inlet transition length, when a stable 
effective thermal conductivity, for the packed temperature profile has been established, then for 
bed and an additional wall resistance. The wall 


; constant heat input is a constant independent 
resistance absorbs over a small radial distance the da 


100 
Ft*Hr 9100 
q 
4580 


| 


LOCAL HEAT TRANSFER COEFFICIENT, h 
5 


- , DISTANCE FROM INLET, inches 


1G. 4. Local heat transfer coefficients along tube length showing effect of flow rate. 
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Fic. 5. Temperature distribution across tube for constant 


k. showing parabolic contour, and temperature drop near 
the wall due to h,.. 


of x. For constant k,, c, G, the equation (2) may 


be integrated to give 


— (3) 


whence 


(3a) 


r*) (4) 


Thus the radial temperature distribution is a 
the 


stream at x will be equal to the fluid t mperature 


paraboloid and the mean temperature of 


at a radius a \ 2. 


(5) 


2k, da 4 
If the heat input is Q B.Th.U. 
from the inlet to the position r, then 


hr. for the area A 


If the total resistance to heat transfer is consi- 
dered to be the sum of the packed bed resistance 
and that of the wall vicinity 


Q 1 
Ah (7) 


(8) 
whence from equations (6), (7), (8) for the mean 


coellicients from the inlet to position Zz: 
(9) 


Thus if the air temperaturse measured at different 
radii immediately downstream from the packing 
outlet are plotted as a function of r*, and the 
best straight line drawn through the points, then 


and also the value 


this construction will give @,, 


of )@ \(r*) for use in equation (8a), Thus with 
values for the overall coeflicient h from equation 
(1) for the whole tube, the effective conductivity 
k, and the wall coeflicient h. ean be obtained 
from equations (Sa), (9). 

The the 


coeflicient h averaged with respect to length was 


arithmetic mean value of overall 
considered appropriate for use in these calcula- 
tions. 


Discussion OF RESULTS 


Experimental data for various flow rates and 
heat inputs have been analysed to give h, h,, 


and k,. As shown in Figs. 6, 7 and 9, each of 
these variables is proportional to G. For turbulent 
flow with Gd » > 100 it is generally recognized 
that the effective diffusivity of the fluid phase 
is proportional to the fluid velocity. This is 
supported by both experimental evidence and 
theoretical analysis [1, 17, 22}. The effective 
thermal diffusivity and hence the effective thermal 
conductivity of the packed volume may therefore 
be expected to follow the same relationship, 
particularly at the higher flow rates when the 
effective conductivity of the solid packing is 
small compared to that of the continuous fluid 
phase. 

The only data for comparison of effective 
conductivity free from the effects of the wall 
resistance are those of CoperLey and MARSHALL 


= 
| | q 
| 
Or 
29 
r 
5 
1956 
4 da/ dr)’ 
Integrating again with respect to 7°, 
0, — 0 
2k, de 2 
: 
G-e 0 a 
Q G-c-« 
A 2 «(a 
( a 
A 4h, 
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© PRESENT DATA TUBE NO.! 
@® PRESENT DATA TUBE NO 2 
~ “VERSHOOR AND SCHUIT TUBE SOmms DIAM 
‘VERSHOOR AND SCHUIT TUBE 43 mms DIAM 


4. 


5 6 7 6 9 10 


REYNOLDS NUMBER 


Fic. 6. 


flow rates. 


(6), Fenix and [81], 8a], BuNNELL et al. 
[4], Scnuter et al. [18], and Brorz [3], shown 
in Fig. 9. The present data are notably lower 
than those of Coperty and Marsa, and very 
much lower NEILL. 


These two groups of workers actually used the 


than those of and 
same apparatus and their results correspond 
closely. The difference shown in Fig. 9 is due to 
the use of the correlating equation proposed by 
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Variation of Nusselt number, based on overall coefficient h, with varying 
(Tube No. 1 and 2 were two constructions of the same dimensions.) 


Fevix and to express their data, which 
equation includes an inverse function of tube 


diameter. 
COBERLY 
values from their experiments. 
of Fevrx and of the considerable influence 
of the tube diameter is difficult to appreciate. 
Brorz [3] on the other hand differs from the 
other authors in that his proposed correlation 


The points representing the data of 


and MARSHALL are however actual 


The suggest ion 
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WALL COEFFICIENT hy 


~+-+-BROTZ 


ux aNO NEILL 


FELIX AND NEILL EQN 


@ PRESENT WORK SPHERES d= O°173in 
© COBERLEY MARSHALL CYUNDERS “@* “ein 
X COBERLEY MARSHALL CYLINDERS 44 * “4in 
SPHERES 2,3, 8mm 


_(CELITE) 


Fic. 7. 


requires the effective conductivity to be pro- 
portional to the particle diameter, which is 
contrary to all other opinions, 

The present data may be represented by the 
equation 


k 0-24 


0.000109 Gd ja. 


If it assumed that the first term represents the 
molecular conduction contribution of both solid 
and fluid, and the second term the convective 
contribution of the fluid, this leads to a value 


for the turbulent Peclet number of. 


c 


0-00109 


9-7 


k 


which corresponds with theoretical estimates and 
experimental values from other sources [1, 2, 17] 
meluding data from mass transfer systems. 

The values of Conerty and for the 
wall coeflicient on Fig. 7 agree with the present 
the 
obtained from a much larger tube. The same is 
true of the data of Feuix and Neri, but the 
correlating equation proposed by these authors 


data for seme mass velocity, although 


predicts values too low due to the inclusion of a 


function of the tube diameter. The scatter of the 


3 


MASS VELOGITY, 2 


Dependence of wall coefficient h,, on mass flow rate. 


points on a logarithmic plot as given by CoperLy 
and MARSHALL is such that their data may equally 
well be correlated by a line of slope unity, corres- 
ponding to h,, o G, as well as by their line of 
slope 4, for h, « G*. Their contention that the 
wall coefficient should be proportional to the 
one-third power of mass velocity, corresponding 
to the normal relationship for streamiine flow, is 
attractive but only tentative, and seems pre- 
mature. The present data cover a much wider 
range of flow rate and indicate the direct depend- 
ence of the wall coeflicient on the fluid velocity. 
It may be noted that the data of Fexrx [8a] and 
PLAUTZ [16] plotted by HANRATTY [10] show a 
mass velocity power of 0-5 for cylinders and 0.77 
for spheres. The equation for the straight line 
drawn through the data on Fig. 7 is 


h 0.04 G 


u 


The present data on Fig. 8 differ from those 
leading to the equation of Hanrarry [10] for 
glass spheres 


hd, 
0-12 


The divergence should not be ascribed to a 
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4, ner 
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63 PRESENT WORK 


} ux 


A 


J 
Fic. 8. Comparison with data of Feirx and PLautz quoted 
by Hanratry [10] for glass spheres. Line A is from 
equation h,, = 0-04 G. Line B is the HANratry equation 


Le} 100 


= 0-12 


dependence on D/d,, as the present packing had 
a value 10-6 whereas the Hanratry correlation 
includes data for D/d, = 24, 16, 12 and 11. 

If it is accepted that the wall coefficients and 
the effective conductivity are both simple linear 


functions of the fluid velocity, the same direct 
relationship with G will be expected for the 
composite overall coefficient h. Though earlier 
data of Leva [14], Co_spurn [7], Cuu and 
Storrow [5] gave a relationship near this, their 
proposed correlations, which provide for a con- 
tinuously decreasing coefficient, h, with decreasing 
G on this account, must be considered inaccurate 
in part and misleading. The coefficient as defined 
undoubtedly does not decrease to zero at zero 
flow rate but to a minimum limiting value. This 
fact was first pointed out by Verscnoor and 
Scuuir [21], and the failure to detect this limiting 
value may be ascribed primarily to errors in the 
measurement of the outlet fluid temperature. The 
practical difficulties are considerable, particularly 
at low flow rates, the experience of Verscnoor and 
Scuvurr being emphatically confirmed in the pre- 
sent work. As the outlet temperature approaches 
equality with the wall temperature, small errors 
of measurement become magnified in the subse- 
quent calculations leading to increasingly low 
values of the coefficient h with decreasing flow 
rate. The data of Verscnoor and Scuutr plotted 
on arithmetic co-ordinates in Fig. 6 are adequately 
represented by straight lines, the intercept and 
slope depending on tube and particle diameters 
and upon the thermal conductivity of the solid 


packing. This simple linearity cannot however 
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Dependence of effective conductivity k, on mass flow rate. 
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be accepted as a precise function, for secondary 
deviations must be « xpected over different ranges 
of flow rate on account of 
(a) the inere asing importance of heat movements 
through the solid packing at low velocities, 
(hb) the change in the nature of the bulk flow from 
streamline to turbulent within the transition 
range, 
changes in the flow pattern on the rear faces 
of the packing particles, 

In normal conditions of operation and when 
Gd 
nature, 
the 


accepte d as a 


100 these effects will only he of secondary 
and a straight line relationship between 


overall coeflicient and flow rate may be 


good approximation. Further 
experi ntal work and a satisfactory theoretical 
an adequate 


analysis are required to prov ict 


correlation of the changes taking place at low 


flow rates 


CONCLUSIONS 


It has been shown that the characteristics of a 


ANDERSON STORROW 


packed tube for heat transfer from the tube wall 
to a fluid moving through the packing may be 
adequately represented by a uniform effective 
conductivity of the packed space and a significant 
wall coeflicient to take account of changes in 
porosity and velocity distribution near to the 
wall. Over the range of air flow rate tested both 
these quantities have been shown to be linear 
functions of the fluid’s mean velocity, and may be 


written for this experimental system as 


10 Gd 
(30 


= 004G 3500) 


k, = 0-24 + 1.09 
Gd yp 1100) 


local cocetlicients 


of the 


showed 


Measurement overall 
of heat that no 
variation occurred with position along the length 
of the tube, 
existed, probably 


transfer important 
but that small secondary effects 
due to inhomogencities of 
packing and velocity distribution, Over a short 
inlet section equal to about six diameters the 
anticipated increase in the overall coetlicient was 


observed. 


NOMENCLATURE 


area of wall for heat transfer ; sq. ft. 

outer radius of packed bed ; ft. 

specific heat of fluid at constant pressure ; 

B.Th.U. Ib. 

diameter of packed bed 2a; ft 

diameter of packing particle ; ft 

mass flow rate of fluid through the packing based 

on the empty tube ; Ib. sq.ft.hr 

overall heat transfer coeflicient ; B.Th.t 
F. 

wall coeflicient for heat transfer; B.Th.| 

thermal conductivity of fluid ; 
hr. F. 

effective thermal conductivity of the packed space; 
B.Th.U.ft. /sq.ft.hr. F. 

total length of the packed bed ; ft 

B.Th.U. hr. 

heat transfer rate for unit length of tube ; 
ft.hr. (B.Th.t 


sq.ft.hr. 


sq.ft hr. 


B.Th.U.ft. sq.ft. 


heat transfer rate ; 
B.Th.U. 
in. hr. where stated). 


q heat transfer rate for unit area of tube wall ; 
B.Th.U. ft.*hr. 

qd heat loss for unit length of tube ; 
(B.Th.t 
linear dimension of packed bed in radial direction ; 
ft. 
linear velocity of fluid based on empty tube ; 
ft. hr. 
total mass flow rate ; Ib. hr. 


B.Th.U., ft. hr. 


in. hr. where stated). 


linear dimension of packed bed measured from 
fluid inlet in the axial direction ; ft. 

temperature of fluid at specific r, 2 position ; F. 
ging at specific circumferential 


temperature of lag 


and a position ; Ik 


temperature of tube wall at specific 2 position ; “F. 


mean temperature of fluid stream ; @,,, at position 


ma 


O ng Ata A 


temperature at wall from extrapolation of radial 


v from inlet to heated section ; 


distribution of temperature in the packing ; °F, 
mean temperature difference : F 
ft? hr. 


thermal diffusivity of the fluid = kc ; 


density of fluid ; Ib. /cu.ft. 
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Abstract—Holding-time distributions are of great importance in a considerable number of 
separation and conversion processes. In the majority of cases it is desired that the distribution 
be as narrow as possible and there are accordingly many examples where a peak signal at the 
entrance gives rise to a comparatively narrow band in the effluent, which has then Gaussian 
shape. 

Several mechanisms have been proposed to account for such a phenomenon. In the present 
study, these mechanisms are analysed on a statistical basis in terms of their basic variables. 
It is also shown that various mechanisms may be combined by adding variances, diffusion 
coefficients or heights corresponding to a theoretical plate. 

It is obviously not possible to identify the mechanism by observing a single Gaussian elution 
curve. Neither is this possible by observing a series of such curves for systems or columns of 
different lengths, in all theories the width of the Gaussian curve increasing in proportion to 
the square root of the length of the system. In the experiments, variation of flow rate, particle 
size, and diffusion coefficients is necessary in order to reach conclusions regarding the mechanism 
or mechanisms responsible for the elution curve. 


Résumé—Les distributions du maintien des liquides sont trés importantes pour un nombre 
considérable de processus de séparation et de conversion. Dans la majorité des cas, il est désiré 
que la distribution soit aussi étroite que possible et il y a plusieurs exemples qui représentent 
une pointe-signal & entrée provoquant une bande relativement étroite dans le courant sortant, 
laquelle bande a la forme d'une courbe de Gauss, 

Plusieurs mécanismes ont été proposés pour expliquer ce phénoméne. La présente étude 
renferme lanalyse de ces mécanismes a base statistique suivant leurs variables fondamentaux. 
On démontre également que différents mécanismes peuvent étre combinés en additionnant les 
variances, les coefficients de diffusion ou les hauteurs relatifs A un plateau théorique. 


Il est évidemment impossible d identifier les mécanismes en examinant une seule courbe 
d'élution de Gauss. I] n'est pas plus possible en examinant une série de courbes similaires pour 
des systémes ou colonnes de longueur variée, toutes les théories impliquant une largeur de la 
courbe de Gauss qui s'accroit proportionnellement a la racine carrée de la longueur du systéme. 
Les épreuves exigent une variation du débit, de la grosseur des particules et des coefficients 
de diffusion pour obtenir des conclusions sur le mécanisme ou les mécanismes donnant la courbe 
d‘élution, 


1. INTRODUCTION publication by Suerwoop [2] for the particwar 


Tue study of times (or volumes) of passage is case of continuous-flow chemical reactors. 

important in a great many processes carried out In such processes it is desired to reduce mixing 

on a large scale in technical units, or on a small in the flow direction of materials that have a 

scale for analytical purposes. concentration or temperature gradient in that 
We refer to Danckwerrs [1] and to a recent direction. 


*Based on an Introduction by the first-named author at the Informal Discussion on “ Mass Transfer in Fluids by 
Diffusion and Convection,” organized by the Chemical Engineering Department of the University of Cambridge 
(July, 1954). 
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This is a general requirement in_ transfer 

processes of the countercurrent and percolator 

typet. Moreover it is important in processes 
involving a single phase, if the composition thereof 

is changing during flow, e.g. as a result of a 

chemical reaction. In the present paper we shall 

not consider countercurrent processes. 

The ideal state nearly always seems to be a 
constant holding time. As regards the transfer 
of matter, the proviso must be made that in 
separation processes it is desired that the holding 
time (volume) be different for the various consti- 
tuents of a mixture. It should, however, still 
be —-as nearly as can be achieved — constant for 
each individual constituent. 

Exceptional cases, where a constant holding 
time is not desired, are the following : 

1. A_ chemical 
stirring, all conditions are constant through- 
out the reactor space. This type of reactor 
designated by Piret as C-star reactor,” 


reactor where, by vigorous 


by others as “ ideal mixer” — is advanta- 
geous when simultaneous reactions are to 
be avoided while consecutive reactions are 
harmless. 

2. Mixing apparatus, designed to decrease in a 
continuous flow the fluctuation of a com- 


position with time (compare Danckwerts 
{1}). 


Examples of processes in which the holding 


time should be constant are : 


1. The heating or cooling of a packed bed by 
means of a Here the 
“ break-through ” is not quite sharp, as a 
consequence of which the thermal efficiency 
This is also 


current of gas. 


of the process is not 100%. 
apparent in the case of heat regenerators. 
2. Various leaching, adsorption ion- 
exchange processes, taking place in fixed 
beds, sharp break-through is 


desirable for the same reasons. 


where a 


3. Some physical separation processes where 
various solutes move at different velocities 
through a column, i.e. chromatographic 


+For nomenclature to be used in characterizing transfer 
processes, see KLINKENBERG [3]. 


Holding-time distributions of the Gaussian type 


processes, based on adsorption, extraction, 
absorption, or ion exchange. In this class 


belongs extraction in  Craiac’s 


in essence this apparatus is a 


also 
apparatus ; 
series of separating funnels [ 4, 5]. 

4. The use of a pipeline for the successive 
transport of batches of different, miscible 
materials, which should collected 
separately. 

5. Chemical reactors, where it is desired 

that all parts of the 


Except 


feed receive the same 
for the 


the reaction will slow 


treatment. rare case of 
zero order reactions, 
down when the reaction partners become 
exhausted. Thus a given reactor volume is 
used less efficiently if some portions of 
the feed stay for too short, and others, in 


consequence, for too long a time. 


The quantities holding time (¢) and holding 
volume (V) are often equivalent in the treatment 
of continuous-flow However, if the 
throughput is not constant it is preferable to use 
the volume, if molecular diffusion or 
heat conduction are to be considered, the time 


pre CESSES, 
whereas, 


is the more useful variable. 

In each case in the following discussion a peak 
of heat or matter may be introduced at 
0, and the distribution is studied 


For other initial conditions the 


load 
Oor V 
in the effluent. 
behaviour of the system may then be easily 
found by carrying out the appropriate summation 
or integration, unless the equations are nonlinear 
(example ; most cases of adsorption chromato- 
graphy). The discussion will be limited to linear 
phenomena. 

The distribution obviously lies between two 
extremes : 

a. All holding times are equal (“ piston flow ”’ 

or “ plug flow” or “ rod-like flow ”’). 
b. The holding time is either zero or infinitely 
high (complete “ channelling ’’). 

In some cases the distribution is described 
by a simple mathematical formula, e.g. in 
the case of continuous flow through a 
single well-stirred vessel. Here, the material 
leaving between ¢ and t + dt after entrance 
into the vessel is the fraction e’ dt of 
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the total amount, where + = mean time of 
passage. 

The three cases described above are shown 
in Fig. 1. 

For such a stirred vessel the distribution curve 
is very w ide. In the above-mentioned processes 
the spread of the distribution is often much 
smaller. It may in fact often be approached by a 
Gaussian probability curve around a mean value. 
It is the purpose of the present paper to analyse 
and compare the various mechanisms that give 


rise to such Gauss curves. 


9 CONVERGENCE TOWARDS rue 
PROBABILITY FUNCTION 


In our further discussion, statistical con- 
siderations will be often used. For an « xplanation 
of basic statistical concepts, like binomial, 
Gaussian, and Poisson distributions, and mean, 
variance, and standard deviation the reader 1s 
referred to the standard textbooks. 

If the distribution in the time of passage is the 
sum of independent individual distributions 


from various successive or simultaneous causes 


FREQUENCY 


and F. SveNITZER 


provided a certain criterion of mathematical 
convergence [6] is applicable —then, as the 
number of individual steps is increased, the 
normal, or Gaussian, or probability distribution 
will always be approached, 

These individual steps may be on different 


scales : 


1. The seale of molecular mean free path. 


Examples : molecular diffusion, 
thermal conductance. 


» <A scale, large with respect to molecular 

mean free path, but small with regard to 

apparatus dimension(s). 

Examples: Eddy diffusion and eddy con- 
ductance in packed beds. 
The same in turbulent flow 
through a long unpacked tube 
(here the scale is small with 
regard to the length). 
Spreading caused by finite- 
ness of transfer rates in packed 


beds. 


HOLDING TIMES 


Fic. 1. Some cases of curves showing the distribution of holding times. 


a. Approaching plug flow. b. Extensive channeling. ¢. One ideal mixer. 
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Holding-time distributions of the Gaussian type 


3. The scale of the apparatus dimension. 
Intermittent 
Craic’s extraction apparatus 
[4. 5]. 


Nonintermittent operation in 


Examples : operation in 


a series of continuous-flow 


mixcrs, 

Sometimes several factors are co-operating 
(e.g. molecular and eddy diffusion). 

In some cases the distribution for the individual 

the 


Examples : 


steps is known and summation may be 
earried out. 
1. Intermittent distribution in stages (CRAIG) : 
fraction p of the solute moves to the next 
stage, fraction 1-p does not move; the 


summation gives binomial distributions. 


2, Nonintermittent operation in stages : the 
individual distribution is according to an 
the 


Poisson distributions. 


exponential law ; summation gives 

In these cases, for a large number of steps, 
the binomial and Poisson distributions may be 
transformed into Gauss distributions. There are, 
however, methods of deriving the Gauss distribu- 
tions directly, based on the general principle 
that the mean equals the sum of the individual 


means. and that the variance | square of the 


standard deviation) equals the sum of the 
individual variances. 

The literature on the subject reveals that 
many mechanisms have been postulated to 


account for observed Gauss distributions. From 
the width of the Gauss distribution one may 
derive a parameter for the particular mechanism 
chosen. Flow in a packed bed may be then 
characterized by a certain equivalent number 
of stages, a number of transfer units, an eddy 
diffusion coefficient, and so on. This, however, 
is no proof that the mechanism chosen is the 
one that in reality caused the observed distribu- 
tion. In view of this state of affairs one must 
examine how the various distributions depend 
on the process variables (time, linear velocity, 
particle size. diffusion coeflicients, ete.), and 
vary such quantities in the course of the experi- 
ments in order to be able to arrive at a conclusion 


regarding the mechanism. Such conclusion can- 


not be obtained from a single experiment only, 
ner can such a conclusion be reached by varying 
the length of the system only, because in the 
region where the Gauss function Is a good approx!- 
mation. the influence of length ts, as we shall see, 
was therefore con- 
the 


mechanisms lead to Gauss functions. 


in all theories the same. It 


sidered useful to summarize how various 
kor such 
comparison one must agree on a standard method 
of expressing the width of the distribution. This 
is chosen to be the ratio of the standard deviation 
o to the mean, p. 

It is also possible that several mechanisms 
that the 


principle mentioned earlier, the squares of the 


co-operate. In case, according to 
ratios due to the individual mechanisms must 
be added. 

Other quantities proportional to the variance 
are also additive. Of particular importance In 
this respect are the (apparent) diffusion coetlicient 
and the height corresponding to an ideal stage. 


The the 


should be expressed in the same units, but it Is 


standard deviation o and mean 


immaterial whether these units are numbers of 
stages, numbers of transfer units, lengths in a bed, 


volumes of effluent or times. The nature of the 


unit will be indicated with a subscript s, tu, a 
v or ft. 

It is recalled that the standard deviation o 
is the the curve at an ordinate 


half-width of 


equalling ¢ 0-607 times the ordinate at the 


mean. Its relation to other quantities which 
characterize the width of a Gauss curve is shown 
in Fig. 2. 


3. Discussion oF InpivipvaL MECHANISMS 


The individual mechanisms now to be discussed 


can be grouped as follows: 


a. Intermittent operation of contacting stages. 

b. Nonintermittent operation of stages. 

ce. Axial diffusion or conduction. 

d. Eddy diffusion and eddy conduction in the 
moving phase in packed beds. 

e. Finite mass or heat transfer coefficients in 


packed beds. 


Furthermore a more distantly related 


mechanism is discussed : 
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POINT OF INFLEXION 
_.» 20.607 


OeSTANDARD DEVIATION 


=O 368 


20 
OV2 180 


Fic. 2. Method to characterize the width of a Gauss curve. 


f. Flow through pipes, when accompanied by vey 


vy vy 
roug (2) 
K 
a. Intermittent operation of contacting stages. 
thus moves on, always and everywhere. The 
A series of equal stages is studied, each con- 
: “ae quantity v is the “ equivalent” or “ effective 
; taining a stationary phase (volume v,,) and a y on 
a : ra volume of the stage. The fraction g = 1 — p 
5 movable phase (volume v,). The stages are “ai 
: remains in the stage. 
numbered in the direction of the flow with a 
. The following distribution is thus obtained 
number 1, 2, ....so that the total 
(Table 1): 
; number of stages, up to and including stage - 
number n, is nm + 1, 
We assume that stage zero initially contains a 
one unit of solute, with a distribution coeflicient Number r of fluid _ 
elements Number of stage n, counting 1 
: i (= (1) or operations from inlet stage = 0 
0 ae 38 4 
; and that the other stages are free of solute. = ae 
. All further portions of movable phase when 0 
as introduced in stage zero are solute-free. 1 P 


1 
2 
It is further assumed that from each stage, 2q 
ry each time it Is used (in each operation ), 4pq® | spy | 
the total amount of movable phase is passed on | 
to the next. Of the solute in the stage, the 


fraction In general the probability P,, of finding a 
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Js (7) 
Ms rp 


While o, itself increases with Vl: 


solute molecule after r operations in the stage 
number n (the (nm + 1)th stage) is: 


> 
ny! (3) 


In operation 1 the mean location of an amount 
q in stage 0 and an amount p in stage 1 is located 
at stage p. Thus, in each operation, for each 
stage, the mean location of the solute is displaced 
by p stages; this is therefore also true for the 
total amount of solute. Thus, after r operations 


the means location of the solute is in stage number 
= rp*. (4) 


The variance is defined as the sum of the squared 
After 
r 


deviations with respect to the mean. 
operation 1, for instance, this is q 

amount q in stage 0, and p x q® for amount p 
in stage 1. Thus, in operation 1 the variance 


2 2 
o, is increased by p*q + q° p 


pq, and after r 
operations the variance of all solute is rpqt. 
By definition, the standard deviation is the 


square root of this variance, thus 


Os, rpq. (5) 


It can be shown that this binomial distribution 
approximates a Gauss distribution — obviously 
with the same p, and o,— when r is sufficiently 
> (p — 4). 


may therefore 


high, viz. when o, 
In that the 
distribution as a continuous function of the stage 


case we write 


number 7: 


where the factor 
f(n)dn = 1. 


We note that the ratio of standard deviation 
to mean location is 


*The subscript s indicates that » and o are expressed 
as a number of stages. 

+The reader may easily verify this, e.g. for r 
working out the sum of quadratic moments ¢ (2p)* 
2pq (1 — + q* (2 — 2p). 


2 by 


ratio 

Hs 

decreases with 4/r. This leads to the familiar 

conclusion that it is advantageous to subdivide 

a given total apparatus volume, as follows : 
The maximal concentration of the solute leaves 

the last of the (n 


of operations r satisfies the condition 


1) stages when the number 


rp bs n-+ 1. (8) 
Then the total volume of effluent is 


n+] 
My vy vy. (9) 
Inside the apparatus the mobile phase travels 
faster (by a factor 1, p) than the solute. Therefore 
the standard deviation in the effluent. expressed 


as a volume, is 


(10) 
agp. 
t p | 
In view of eq. (8) we obtain 
o, = 1 4/q(n +1). (11) 
p 
Thus the ratio 
q (12) 


In those cases where g +1 experimental 


vield immediately the number 
1 (Mayer and Tompkins (7}). 


values of o, 
of stages n 


see also our equation (18). 


b. Nonintermittent operation of stages 
Equations for nonintermittent operation of 
stages are obtained by letting p +0; ¢g +1; 
r ow, wherees rp remains constant. 
Since p is the ratio of the element by volume 
which is displaced to the equivalent stage volume, 
the quantity 


x rp (13) 


represents the total volume passed, expressed 
in the equivalent stage volume as the unit. 
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Since n «< r we may transform eq. (3) to 


(14) 

This is Poisson’s distribution law, which ts 

well known for systems consisting of one phase 

(Ham and Cor [8]; and WeBer [9]; 
Jounson and Epwarps [10]). 

The derivation of the Gauss Law as given 

ander (a) is modified accordingly by replacing 


q by 1 and rp by 2 so that 


(15) 
o? (16) 
o 
(17) 
Hs Ni 
while at the end of a column consisting of n +1 


st: oes, 
We may compare the variance of the binomial 


distribution and the variance of the Poisson 


distribution, for an equal number of stages, 
with the same amounts of the two phases in the 
stages and with the same total throughput. 

In this case the mean, rp, for both is the same 


and the ratio of the variances 1s: 


binomial q 3 (18) 
oO” Poisson 


It is interesting to note that this is the same 
factor, suit: bly tri nsposed, which GLUECKAUI 
[11 obtained when the effects of the two types 
of stages wer compared by another method, 
Here the apparent diffusion coefficients, required 
to account for the presence of the stages in 
equations for diffusion, were in the same ratio 
to one another. 

Since the factor q depends on the kind of solute. 
for an actual separation the two methods of 
calculation give essentially different results 
(GLUECKAUF (11)). 

The statistical treatment may be extended to 
cover the case where the solute is introduced 
over a period of time. The Gaussian curve will 
be approached if the variance in the feed is small 
enough with regard to the variance produced by 


*P, . is a maximum for 2 n. 


and F. 


the stages. A mathematical expression for this 
has been obtained by Vax Deemwrer, ZUIDERWEG, 
and KLINKENBERG [12] who also show how the 
standard deviation may be derived from concen- 


tration curves obtained with a longer feed. 


c. Arial diffusion or conduction 


At any given instant the distribution as a 
function of height is given by a Gauss function, 
The width of this function is derived as follows: 

In a column of height 1, through which fluid 
is passed with an interstitial velocity u, the 
nominal time during which a solute diffuses ts 
the time of passage 7 of that fluid, even 

u 
if the solute is retarded by a stagnant phase. 

The effective diffusivity, which is diffusivity 
corrected for tortuosity, will be indicated by 
D. According to well-known law, 
see e.g. GLASSTONE (13), the variance 1s 2Dr. 


so that: 


» 
of (19) 
and 
pe? (20) 
ul 


Although at a given instant a graph of concen- 
tration vs. height ts symmetrical, the graph of 
effluent concentration vs. time will be skew 
because the various clements of solute have been 
diffusing for different times. If o, l the 
skewness may be neglected. 

For conduction of heat analogous equations 
are easily derived by replacing the molecular 
diffusivity D by the thermal diffusivity a. 

Here an effect must be considered which, 
in the case of mass transfer, seems to be of 
subordinate importance only, viz. transport 
through the particles. 

Bakuvurov and Boreskoy [1 have examined 
apparent diffusivities for heat and matter in 
packed beds through which a gas ts flowing. 
For transfer of heat and mass in a bed of crushed 
porcelain the two diffusivities were equal, while 
in a bed of iron shot heat moved more quickly 
than matter. This difference is almost certainly 
due to the higher thermal conductivity of iron. 

In these experiments eddy conduction and eddy 
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diffusion, to be discussed sub d, played a role ; 
however, these mechanisms must have had the 
same effect in both cases. 

Method of 
conductiv ity of packed beds have been examined 
by Sincer and [15]. 


computing an overall thermal 


d. Eddy diffusion and eddy conduction in the 
moving phase in packed beds. 

In some cases the flow in a porous bed is so 

1) that 

the inertia terms in the basic hydrodynamic 


slow (Re, calculated on particle size 


equations (the only terms which are not of first 
order in the velocity ) may be neglected. 

1 the pattern of the 
multiply 


In the domain Re, 
streamlines remains the same if we 
all linear dimensions or all velocities by a constant 
factor. 

If the solute stays with the liquid in which it is 
contained — i.e. if there is no diffusion, and no 
transfer of mass — the following reasoning shows 
that flow through a porous bed causes a spread 
in residence time which obeys the laws of diffusion. 
The mechanism is therefore called eddy diffusion 
by analogy with the well-known eddy diffusion 
in turbulent flow. 

In view of the similarity of the flow pattern 


for Re 


irrespective of the velocity u (interstitial velocity ) 


in geometrically similar packings, 
and the particle size d,, upon passage through 
one layer of particles (thickness d,) the standard 
deviation of the time of passage must be a constant 

P itself. Or, 
u 

in terms of variance, we may put : 


fraction of the time of passage 


of = 2A (21) 


i.e. for the total packed bed with approximately 


layers : 
Pp 
2Ald 
a (22 
u? 
Furthermore 
l 
u 


*The reason for using the factor 2 is seen later. 


so that 


d,, (28) 
l 

By analogy with eq. (20), since the left-hand 
members are the same, it is seen that the quantity 


E A d, u (24) 


now takes the place of the molecular diffusivity 

D. This quantity is called the eddy diffusivity E. 

Combining equations (23), (24), and (17) gives 

a relation between the eddy diffusivity E and 

the equivalent number of theoretical stages of 
the packed bed, n +1: 

ul 
E (25) 


2(n + 1) 


This result derived independently by 
Reman and Brurnzee. [16] and by Kramers 
and ALBERDA [17]. It is also closely related to 
GLUECKAUF’s treatment [11]. 


In eq. (23) it is seen that the velocity of flow 


was 


influence on : 


has no all that is important 


is the ratio of length to particle diameter d,, 
together with the constant A which characterizes 
the type of packing. 
Combining equations (24) and (23) gives a 
between the equivalent number of 
1 and the ratio of height of the bed 
to particle size d 


relation 
stages n 
P 


1 /l 
= 


p 


(26) 
Relatively little has been published about 
values of A for axial eddy diffusion, presumably 
because packing irregularities causing channelling 
have a far greater effect than the flow pattern 
in a carefully prepared packing. 

A. KLINKENBERG in 1945, in an unpublished 
study of the eddy diffusion in the absence of 
mass transfer (water followed by aqueous solu- 
tion, or vice versa, through carefully packed 
quartz sand), obtained A values ranging from 
2 to 4. 

In one typical experiment the sand, a 70-120 
mesh fraction (d, = 125-210 4), was contained 
in a 83cm diameter tube, filled height 55 cm. 
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free volume F, = 38°8%, and the displacement 
was of water by 2-molar nitric acid solution 
(Experiment 1), and vice versa (Experiment 2), 
both flows in downward direction under conditions 
where Re Re 

The effluent-concentration curves the 


following characteristics (Table 2) : 


Table 2. 


Experiment 50°, point Standard 
No. Hy in ml deviation A 


in ml 


1 34 
2 141°8 


t 


Equating the ratio of standard deviation to 
mean in terms of volume to the one in terms 
of time (eq. 23), with d, = 168 10 °m and 
l 55m one finds the A values indicated in 
the last column. 

Judging from the relatively close  corres- 
pondence of the A values, there cannot have 
been much interference by buoyaneys effects. 

KraMers and ALBERDA [17] determined 
\}~1-1 for Raschig rings in a_ packed bed, 
Danckwerts [1] in similar experiments found 
A ~ 18. 


e. Finite mass or heat transfer coefficients in 
packed beds. 

In the treatment of transient beat or mass 
transfer in packed beds following a sudden 
change of concentration at the inlet (say from 
0 to ¢,), it is usual to introduce dimensionless 
co-ordinates Y and Z with the following signi- 
ficance for mass transfer (see e.g. KLINKENBERG 


[18]): 


al 
Y height expressed in transfer units re 
uf, 
Z time since the fluid introduced at the 


moment of the jump reached the point 
under consideration, expressed in time of 


ak KF, 
passage per transfer unit zm § —. 
xl Pi, Fy 
uf, 


*For heat transfer, « is to be replaced by 
pf 


and F. Ssenrrzer 


This Z may thus be seen as the difference of a 
ak KF, 

true dimensionless time and 
Fy Py 

times the dimensionless place co-ordinate Y. 
As a first approximation, for large Y and Z, 
the concentrations in solid and fluid are then 


both proportional to 
VY)] 
The effect of a peak load is found by differentiation 
with regard to Z: 
2 \ 

This distribution assumes appreciable values 
only for small (4/Z — +/ ¥), hence for large 
Y only in the vicinity of Z = Y, while the inde- 
pendent variable 4/Z — 4/ Y may be replaced by 
Z. 

At the maximum : 

\ Z \ \ Hol Z) VY 0. 
At a distance o,,(Z) away from the maximum : 
\ Z \ 4 V bm (Z) on, (Z) \ Y 4 Vv 2. 

Hence the difference of these two Z values is 
in first approximation. 

On (Z) = 2¥ (28) 

Thus, the distribution of Z is characterized by a 
(Z) = Y and a variance of (Z) = 2Y. 


py, 


In view of Z = 7 -!¥Y the distribution 
u 
of T is characterized by a mean y,, (T) 


KF 
Fy, 
Hy, (Z) Y this becomes : 


(T) (1 


Z ). Since Y is a constant and 
y. 

Fy 

On the other hand the variances of the T and Z 


distributions are the same, since 7' —Z is a 


constant, so that 
(T) = (Z) = 


Thus the ratio of o and yp becomes : 


Mu (T) Fy + 
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In view of the physical significance of Y it is 
seen that, again, o increases as the root of the 


height and “4 is inversely proportional to that 


root. 

Kinetic Theories of Mass Transfer. In kinetic 
theories rate mechanisms are proposed for the 
capture of molecules of the solute by the stagnant 
phase, and for the reverse process of release of 
Such 


describe 


these molecules. theories have been 


developed to adsorption and ion 
exchange. 
If both of the 


equilibrium isotherm is linear and the theory leads 


reactions are first order the 
to the same results as the one discussed sub e. 
This is a special case of the theory developed by 
THOMAS [19] and others, see also GippINGs and 
Eyrine [20] and Kunkenpere [20a]. 
Reactions of second order and 
first 
by Tuomas [19], Ginuitanp and Bappour [21], 
Bappour, D. Gotpsrern and Epstein [22], 
S. Gotpstein [23] and Vermeucen and Hiesrer 


[24] lead to a These 


and second order combined, as examined 


non-linear isotherm. 


theories are consequently outside the scope of 


the present paper. 


f. Flow through pipes, when accompanied by 
radial diffusion. 


Flow through long pipes in either the laminar 
or the turbulent region is more complicated 
than the cases discussed so far in the sense that 
the concentration in the mixing zone of two 
miscible fluids or in a chemically reacting mixture 
depends on two co-ordinates. 

However, the interplay of differences in holding 
time of liquid or carrier gas with radial diffusion 
or eddy diffusion is such that in long pipes the 
holding time of the solute, averaged over a cross 
section, is again distributed according to the 
probability laws. 

Laminar Flow. Laminar flow has been studied 
by Bosworrn [25], who was interested in the 
influence of the holding time distribution on the 
efficiency of a tubular reactor for reactions of 
gases, and by Tay.or [26]. who studied the 
axial spreading of an injection in fluids. See also 
Van Deemrer et al. [27]. 


reactions of 


Boswortu [25] found that molecular diffusion 
affects the fluid 
moving in laminar streamlines through a cylinder 


residence times of elements 


in two ways, viz. 

1. Radial effects by which the distribution of 
residence times is made more homogeneous. 
Longitudinal effects by which this distribu- 
tion is made less homogeneous. 

the radial effect 
diameter d > 25 


that 
tube’s 


showed 
the 
The case that this effect is not negligible 


Boswortu 
negligible when 
DL u. 
was discussed by Tay.or [26] who combined 
molecular diffusion with the variation of liquid 
velocity in the tube. This resulted in an effective 
axial diffusivity 

d? u* (29) 
192D ~ 

We have mentioned in Part 2 that variances 


due to statistically independent mechanisms 
must be added and shall apply this principle 
in Part 4. 


parabolic velocity distribution and the effect 


The axial spreading caused by the 


thereon of radial molecular diffusion are, however, 
interdependent. 

Turbulent Flow. Turbulent flow 
studied by Tayior [28]. who found that 
particles, released in a turbulent constant flow 


has been 


through a long straight pipe, will spread in axial 
directions from the initial position in the streaming 
fluid according to the laws of statistics. An effec- 
tive axial diffusivity A* was introduced and 
related to the friction factor in the pipe. From 


TAYLOR’s equations we can derive the ratio 
o 
2k * 


ul 
in complete analogy with eq. (20). 


COMBINATIONS OF MECHANISMS 


The mechanisms,giving rise to Gaussian distri- 
butions in the effluent, if occurring successively 
and independently, may be combined according 
to the rules known from statistics, i.e. by adding 
the variances. 

Quantities H.E.T.P. and the 


such as the 
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(apparent) diffusion coefficient, which are pro- 
portional to variance, also have the additive 
property. 

It is not quite clear whether it is permissible 
to add variances if the mechanisms take place 
simultaneously, since it is then more difficult 
to judge whether they are statistically indepen- 
dent. One such combination was studied by 
Vax Deemwrer, Zurmerwec, and KLINKENBERG 
[12]. For this combination (diffusion and finite 
mass transfer rate) it is reported that this is 
indeed SO. 

In order to aid this type of work we have 


prepared the following table (Table 3) : 


Table 3. 


o* 2D 


Fy 2uF;, 
-KF, 


al 


(30) 


If it is desired to derive of (i.e. variance in 
terms of times), the members of this equation 
must be multiplied with the square of the mean 
time + required for the passage of a band through 
the column : 

y= + (31) 
uKF, 


so that 


32 


This is identical with the result derived in [12]. 


Survey of characteristics of holding time distributions. 


1 2 
Mechanism Unit 
a. intermittent stages —_ stage 
b. nonintermittent stages stage 
c. uniform flow and molecular diffusion length 
d. uniform flow and eddy diffusion time 


transfer 


e. uniform flow and finite rate of transfer nit 
uni 


NOTE. 2 


3 5 


Measure for 


Mean Variance relative width 
m o* of band 


q 
rp rpq 
rp 


al (1 al Fyy 2uF, 
uF, Fy uF, Fy t KF, al 


The ratio ae (square of quantities in column 5) is additive. 


By way of example we discuss the combination 
of molecular diffusion and finite rate of mass 


transfer, mentioned above. 
_ is found by adding 
the squares of the expressions in the 5th column 


The relative variance 


of the table, viz. 


If it is desired to express a height of a theoretical 
plate in terms of its contributing factors, the 
survey as given in the table (column 5) indicates 

that such height - is proportional to a * and that 


D 


u 


its contributing terms are proportional to 
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Holding-time distribution of the Gaussian type 


(for molecular diffusion), independent of u (for ¢,; = concentration in mobile phase suitable unit 

eddy diffusion) and proportional to u (for finite = 

= same in stationary phase suitable unit 
rate of mass transfer, assuming « to be indepen- “ 

dent of u). This relationship is examined in a Cy = initial concentration suitable unit 
quantitative way in the publication [12]. per m® 

As a further example we may discuss SAUNDERS D = effective longitudinal (molecular) 
and Forp’s study of the heat transfer by a diffusivity "eee 
current of gas to a bed of solid spheres [29]. agi 7 
d, particle diameter m 
he results were evaluated in terms of a transfer BE longitudinal eddy diffusivity m* ‘sec 
coeflicient (see Part 3-e) and the quantity } Fy, fraction by volume occupied by 
mobile phase 
was found to be propor tional to d (lc. page Fy, fraction by volume occupied by 
305 P). so that. in our symbols, P immobile phase 
hl h heat transfer coefficient J °C"! 
Y = 0-63 (33) ‘ 
pe, u F, u F, d, K distribution coefficient ( in 
equilibrium) 
Thus the heat transfer coeflicient -opor 
us the heat transier coetlicient 2 Is propor- K* = effective longitudinal diffusivity 
tional to the velocity u. in tube m® ‘sec 
When working on the basis of eddy conduction, 1 = length of packing m 
however, since AF, Fy. according to the n number of theoretical stages 
probability 
xl 
“ Survey,” column 5, we must replace BF by p = fraction of solute present in 
un mobile phase 
a\d ° Consequently on this basis SAUNDERS and q fraction of solute present in 
resi oh ‘rpreted as - 
ol sm ts might be rpre Re, Reynolds’ number based on 
stration of the validity of the afore-mentioned particle size 
theory for eddy conduction with a constant value r number of stages 
of A. viz. sa S = total volume of mobile phase 
A O79 flowed through 
limensionle time, ‘fined ji 
his A-value is even lower than the ones given - : = ees time, defined in 
art 5-e 
in Part 3d, as observed for liquids, but on account t time 
sec 
of the argument given earlier in Part 3d, this u = interstitial velocity m ‘sec 
might be expected for gases where there is con- V effective total volume m*® 
siderable interaction by radial conduction. effective stage volume 
Moreover the value of the Reynolds group, vy volume of mobile phase in one 
theoretical stage m® 
calculated on particle diameter, was the 
volume of immobile phase in one 
ange of 10-1,000 whereas the considerations theoretical stage m3 
of Part 3d were deduced for Re 1. a Sv 
It is beyond the scope of the present article : : al 
to examine the details of the mechanism. The uF, 
above considerations, however, may suflice to Z = corrected dimensionless time, 
show that mechanisms other than resistance to defined in Part 3¢ 
mass transfer coefficient per 
heat transfer at the particle surface have to be I 
unit volume of packing sec”! 
considered. A coefficient defined in eq. 21 
List OF SYMBOLS and 24 
Unit m mean location depends on i 
. 
a = thermal diffusivity m/sec subscript 
c concentration (general) suitable unit o standard deviation or context 
per m® p = density kg m3 
c specific heat Jie* C= Tt mean time of passage sec 
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Subscripts with » and a indicate the unit in which they tu = as a number of transfer units 
are expressed. v = as an effluent volume 

1 = as a length Subscripts indicating the phase : 

s as a number of stages I = in mobile phase 

i= as a time Il = in stationary phase 
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Longitudinal diffusion and resistance to mass transfer as causes of 
nonideality in chromatography 


J. J. van Deemwrer*, F. J. Zuriperwec* and A. KLINKENBERG** 
(Received 1 February, 1956) 


Abstract—The mechanisms of band broadening in linear, nonideal chromatography are 
examined, 

A development is presented of a rate theory for this process, wherein nonideality is caused by: 

axial molecular diffusion ; 

axial eddy diffusion ; 

finiteness of transfer coeflicient. 

The correspondence with the plate theory is given, so that the results can also be expressed 
in heights equivalent to a theoretical plate. 

The plate theory has been extended to the case of a finite volume of feed ; the requirement for 
this feed volume to be negligible has been examined and a method is presented for evaluating 
concentration profiles obtained with a larger volume of feed. ‘< 

An analysis is given of experimental results, whereby the relative contributions to band 
broadening for various cooperating mechanisms could be ascertained. 


Résumé— Une étude a été faite des mécanismes de lélargissement des bandes dans la chromato- 
graphie linéaire non-idéale. 

On a développé une théorie sur la vitesse de ce processus, selon laquelle la non-idéalité est 
causée par : 

la diffusion moléculaire axiale ; 

la diffusion turbulente axiale ; 

le fait que le coefficient de transfert a une valeur finie. 

La correspondance établie entre cette théorie et celle des plateaux permet d’exprimer les 
résultats obtenus aussi en termes de la hauteur de colonne équivalente 4 un étage théorique. 

La théorie des plateaux a été rendue applicable sur le cas d'un échantillon de volume fini ; 
les conditions dans lesquelles le volume de cet échantillon est négligeable ont été étudiées et on 
présente une méthode pour interpréter les profiles de concentration obtenus avec des échantillons 
plus grands. 

Une analyse des résultats expérimentaux a permis de déterminer les contributions relatives 
des différents mécanismes en question a lélargissement des bandes. 


1 INTRODUCTION 1. The equilibrium concentrations in the two 


For the design of efficient chromatographic phases are proportional (“ linear 


columns it is desirable to have a background of chromatography ”). 


The exch; ESS 1S > ies , 
theory, in order to know how various factors 2. The exchange process is thermodynamically 
influence the sharpness of separation. reversible, i.e. the equilibrium between 


In the theories of chromatography simplifying particle and fluid is immediate — the mass 
assumptions are introduced in order to make transfer coefficient is infinitely high —- and 
the basic differential equations accessible to longitudinal diffusion and other processes 
mathematical treatment. As the nature of the having a similar effect can be ignored. 
assumption depends on the type of chromato- Such a process is commonly called “ ideal 
graphy adopted, a survey of the possible assump- chromatography.” 


tions should first be made. These are : By making these assumptions or not in the 


*Koninklijke /Shell-Laboratorium, Amsterdam. (N.V. De Bataafsche Petroleum Maatschappij). 
**N. V. De Bataafsche Petroleum Maatschappij, The Hague. 
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two cases respectively one arrives at the following 


four possibilities : 


Ideal Nonideal 
chromatography chromatography 


Linear isotherm I Ill 


Nonlinear isotherm Il IV 


Case I covers the simplest theory still possess- 
ing the essential features of chromatography. The 
retardation of solute with regard to solvent 
depends on the product of the distribution 
(partition) coeflicient between the phases and the 
ratio of the amounts of those phases present in 
the column, but not on concentrations. The 
shape of a band, during its movement, thus 
remains unchanged. Furthermore, different 
solutes behave independently [1]. 

If a band of mixed solutes is introduced in a 
column, it is a matter of simple arithmetic to 
find the requirement for the individual bands to 
become separated ; this point need not be gone 
into further. 

Case II is of particular importance for 
adsorption chromatography in liquids where the 
effect of nonlinearity of the isotherm can usually 
not be neglected, and where it may be assumed, 
on the other hand, that mass transfer is sufliciently 
fast because of the small size of the adsorbent 
particles, and that, in spite of the low linear 
velocity of the fluid, axial diffusion is still 
relatively unimportant. The first author to 
provide an adequate treatment of the single- 
solute case was De [2]. 

The bands change during their passage through 
the column and usually dev elop a Sharp front and 
a long tail. 

A special complication 1s introduced by the 
fact that the solutes may affect each others’ 
movement, so that the case of two or more solutes 
usually cannot be derived by superposition of 
single-solute cases. Strictly spec king, it is not 
even accessible to rigorous treatment [3, 4, 5]. 

Case II is of particular importance for partition 
chromatography with a moving liquid or gas. 
Here the assumption of the linear isotherm is 
usually a good approximation. 
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Case IV, the most general case, is probably 
required for adsorption chromatography of 
gaseous mixtures. Furthermore, certain kinetic 
theories mentioned by KLINKENBERG and 
SJENITZER [6] covering adsorption and ion 
exchange, belong in this category. In_ these 
cases, the mathematics become very involved. 

In this paper a discussion will be presented of 
Class III of the chromatographic theories. 

There are two basically different approaches, 
known from many other transfer processes. In 
the plate theory, the column is conceived as 
consisting of a number of stages or plates in 
each of which there is equilibrium between the 
two phases. The second approach, which may be 
called the “rate theory,” considers the actual 
continuous column. Diffusional phenomena etc. 
are taken into account for describing the un- 


sharpness of a band. 


The Plate Theory 

In the plate theory (7, 8 | starting from the 
assumption of a linear distribution isotherm, the 
separating efliciency of a chromatographic column 
is characterized by the height equivalent to a 
theoretical plate (H.E.T.P.). This H.E.T.P. is 
an empirical quantity and the theory does not 
deal with the mechanisms which determine it. 
It is, however, of much practical value for 
calculating the effects of the column length and 
the method of feeding the column on the separa- 
tion obtained. The H.E.T.P. depends on the 
flow rate and on the nature of the substances 
treated just as it docs in the packed columns 


employed in other transfer processes. 


The Rate Theory 

The rate theory in principle provides all 
information on the influence of kinetic phenomena 
such as rate of mass transfer between phases, 
rate of adsorption or chemical reaction, longi- 
tudinal diffusion and flow behaviour on the 
history of a band in the column. 

Here, the basic difficulty lies in arriving at 
an adequate physical concept of the phenomena 
determining the movement of a solute in and 
around the particles. This is where various 


authors differ in their way of treatment. Such 
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a concept is of necessity an idealization, and its 
acceptability will have to be checked experi- 
mentally. 

One method of treatment, analogous to well- 
known work on heat transfer (see review by 
KLINKENBERG [9]) introduces solely a transfer 
coefficient (e.g. and Boyp [10]). It 
is, however, easy to see that this method lacks 
detailed physical background, since the usual 
for the er two-film 
concept do not apply. 

and Boyp [10] and Boyp, Apamson, 
and Myers [11] have interpreted this coefficient 
by a diffusional mechanism through a liquid film. 


conditions film concept 


Band broadening by a finite mass transfer 
coeflicient bas been combined with band broaden- 
ing the of 
Laripus and Amunpson [12]. 

The use of an empirical mass transfer coefficient 


by axial diffusion in treatment 


has been avoided by Tunrrsku [13], who took 
eddy diffusion and intraparticle diffusion into 
consideration, and by GLUECKAUF [8] who 
related the H.E.T.P. to particle size, particle 
diffusion and diffusion through the film surround- 
ing the particles. GLUECKAUF also showed how 
his theory may be used to explain many experi- 
mental effects in ion-exchange chromatography. 

The present paper may be seen as an extension 
of the work of GLUEcCKAUF and partly as an 
alternative, namely the 
between H.E.T.P. and contributing mechanisms 


in so far as relation 
is observed, 

It also uses AMUNDSON and Lapipus’s approach, 
but with a considerable simplification. It will 
namely be shown that under certain conditions, 
which are of practical importance, AmMuUNDSON 
and Lapipus’s formulae may be simplified to a 
Gaussian distribution function which is directly 
comparable with a 
the plate theory. 
lines of thought leads to a simple relation between 
the H.E.T.P. and various parameters of the rate 
theory. The result will be applied to two examples, 
namely to the 
Simpson and WuHEeaton [14] and to results for 
gas-liquid-partition chromatography obtained by 
Kevutemans [15]. 

Reference is also made to the publication by 


A combination of these two 


ion-exclusion experiments 
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corresponding formula of 


of 


KLINKENBERG and Ssentrzer [6] the subject 
of which is the distribution of retention times 
in flow systems, particularly of the Gaussian 
type. 
2. Tue Tuerory or THE Piate 
Derivation 

In its most simple form the plate theory has 
been developed and applied by Martin and 
SyncE [7]. 


will be extended somewhat 


In this section the plate theory 
to include the case 
where the influence of the feed volume in elution 
chromatography may not be neglected. 

Consider the history and the elution curve of 


a single band in partition chromatography. 


Let v, and v,, be the volumes of the moving and 


the immobile phase in one theoretical plate 


(Fig. 1). In the course of the development of 
as as 

“Zp 
STAGE Nr nel 
Fic. 1. Partition chromatography in stages. 


the theory these volumes will be taken as constant, 
i.e. the concentration of solute in both phases has 
to be low. The effective plate volume v may then 
be defined as 


v vy T Uy] K. (1) 


is assumed to be 


factor K 
independent of the concentrations (linear isotherm ) 


The distribution 


and is given by 


Ke;;. (2) 


The concentrations c,; and ¢,,; are measured in 
suitable units per unit of volume. 

Let the feed have a concentration ¢c, and the 
feed volume be equal to 4. The material balance 
for the first plate, if a volume dS of the moving 
phase flows through it, then becomes : 


vy de; T Uy dey, Cr dS 


Co dS for 0: 
(3) 
0 for S > A. 
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For the other plates 


dey + V de + 4S 
dS (n 1). (4) 


If we introduce (1) and (2), and we write 


y, and w. (5) 
Co v 


Cy, n 


(3) and (4) are reduced to 


A 

j 1 for 0 < w: 
at 7 
‘ (6) 

di 
Yn = Yay 1). (7) 
dw 
The initial condition is 

UP Yo = = 0. (8) 


As may easily be verified, the solution of (6), 
(7), and (8) for all stages including the first is : 


{ 
(nm )! 1 
0 
y * dw' forw >; (10) 
(n 1)! 1 
u 
or 
s 
g’\" 1 
vJ (n—1)! 
for 0<S <A (11) 
s 
Che Co S’/t | ds’ 
v J (nm 1)! v 
S-A 


for S > A. (12) 


These formulae are sufficient for calculating 
the position and shape of a moving band. In 
the following, formula (12) in particular will 
be considered more closely, since it represents 
the situation at the end of the column. 

The integrand of (12) is the Poisson distribution 
function. This function may be approximated 
by a Gaussian distribution function if the argu- 


27 


ment S’/v is large (say S’/v > 100), This will 
always be the case when the number of plates 
is not too small and the band has arrived at the 


end of the column. 
al v 

1 


{ (S'/v —n)*) 


13 


exp 

With sufliciently large values of 8S", v no serious 
error is introduced by replacing the Gaussian 
distribution (13) by 


2n | 


With (14) and by putting x 1 = n equation 
(12) becomes 
S’ /v n)* 
2n | 
S-A 


Cr ds’. (15) 


The elution curve of the band according to 
(15) is indicated in Fig. 2. The concentration is 


a maximum for 


Gan 
- a 
10+ 
POINT 
be 4 Sf - 
ove 
Fic. 2. Elution curve of a band. 
nv + 4A. (16) 


max 
According to the law of conservation of mass, 
the area under the elution curve is equal to the 
feed volume A. 
For small feed volumes (4 < 8S) the band 
obtains a Gaussian shape : 
n)*) 


exp 


| 2n (17) 


Chin 


As will now be shown the condition A < 8 
may be brought into a quantitative form. 
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We found it convenient to define the width and 
AS of the elution curve as the distance between 


A a + 28 
the points of intersection of the tangents in the ; ats ‘ 
inflection points with the horizontal axis (Fig. 2). \/2n 5 (erf — (21) 
F For the Gaussian curve this value of AS is just = = 
- equal to 4 times the standard deviation. where 6 is given by 
2 In order to relate width and height of the bet — (a 4 8) Hata, (22) 


elution curve to the number of plates it is con- 


venient to introduce the dimensionless quantities ’ 
It is seen that height and width 


Ss 4 Co / max 
$= anda (18) AS 
van vyn As = - , are only dependent on the parameter 
Expression (15) then becomes A 
a These relations are graphically 
vyn 
eo dy’. (19) represented in Figs. 3 and 4. 
V 27, Fig. 4 shows that for 
sa 
A straightforward calculation shows that a r. < 0-5 (23) 
Che vyn 
Co} max 24/2 the width AS becomes almost independent of 


MAX. 
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Fic. 3. Height of elution curve. 
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i 
2 3 


Fic. 4. Width of elution curve 


the feed volume. If in this case we denote the 


width as AS,, we have 
AS, (24) 


This result means that the elution curve is 
Gaussian. 

It may further be derived from (21) and (22) 
that 

AS A + vaf/2an for - 3. (25) 

Another interesting characteristic of the elution 
curve may be derived from the above equations. 
Let the retention volume (S,) be defined as the 
volume of mobile phase that has flowed through 
the column from the instant the feed is stopped 
until the moment the solute appears in the 
effluent. The value of S, at which all the solute 
is just stripped from the column is then called 
the final retention volume, S,, (see Fig. 2). 
According to this definition 


ni A T AS. 


If we substitute AS by using eq. (24) and (25) 
it then follows for a > 3 that 


= nv + 0-81 AS). (26) 


Since AS, is constant, (26) shows that the final 
retention volume is independent of the feed 
volume. The final retention volume and AS, 
being easily obtainable experimentally, eq. (26) 
is very useful for calculating the total stage 
volume nv. The subsequent calculation of the 
number of plates from the width and height 
of the elution curve is then most easily performed 
with the aid of Figs. 3 and 4. 


Number of plates needed for chromatographic 
separations. 

The chromatographic technique is generally 
used for separation of two or more solutes. There- 
fore the question arises how large the minimum 
number of plates is for obtaining such a separation. 
Obviously the degree of separation depends on 
the overlap of the elution curves of the solutes 
emerging successively from the column. Since this 
overlap decreases if the elution curves of the 
solutes become narrower, it immediately follows 
that for a given degree of separation the minimum 
plate number is obtained when small feed volumes 
(a < 0-5) are used. In order to obtain a quantita- 
tive expression for the minimum plate number, the 
degree of actual separation that results has to be 
defined. A general solution to this problem 
has been offered by GLUECKAUF (16). However, 
a simple case is obtained when the tangents of 
the elution curves of the subsequent solutes 
just touch at the base of the chromatogram. 
With v, and v, being the effective plate volumes 
of solutes 1 and 2, the distance between the 


maxima of the bands is 


— NV, 
or 
+1\? 
(27) 
vw, 
when ’ (28) 
* 


» is the ratio of the distribution coefficients ; 
in vapour-liquid operations it is called the relative 
volatility of the solutes. Eq. (26) is valid only 
for the specified degree of separation. A calcula- 
tion shows that, for the case of an equimolar 
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of the eluted 
» if a separation is effected 


mixture of solutes, the purity 


material is about 9 

at the minimum between the two concentration 

peaks, 
It 


plate number needed for chromatographic separa- 


is interesting to compare the minimum 


tions with the minimum of plates needed in a 
for batch 
distillation. Fig. 5 shows this comparison for the 


countercurrent — process, instance 


MIN/MUM NUMBER OF 


THEORETICAL PLATES 
7 
1/000 
D CHROMATOGRAPHY 
r 
ATCH DISTILLATION 
42 6 8 é 
REL AT ry 
Fic. 5. Comparison of plates needed in chromatography 


and batch distillation. Separation of equimolar mixtures 


into two fractions of 9 


0 
oe 


» purity. 


separation of equimolar mixtures; the plate 
numbers required are plotted as a function of 
the value of y. It is seen that in chromatography 
much higher plate numbers are required than in 
This 


undoubtedly is caused by the fact that in chroma- 


distillation, especially at low values of ». 
tography only a small fraction of the total plates 
is in use at a given instant, while in distillation 


all the plates are used continuously throughout 


the process. Though it thus turns out that the 
countercurrent process is much more economical 
as regards the number of plates required, it should 
be kept in mind that the high number of plates 
needed for chromatography is much more easily 
realized than the lower number of plates in 
distillation. 


Application to gas-liquid partition chromatography 

The predictions of the effect of feed volume on 
the width and height of the elution curve have 
been checked experimentally. Use was made of a 
chromatographic column 1-80 m long and 6 mm 


in diameter. It was packed with 50-200 mesh 


Table 1. 
Experimental data on width and height of elution 


curve for gas-liquid partition chromatography. 


A A S. (Chm 
ml butane mil ml Co | max 
in feed 
exp. calc.* 
32 10 163 1100 0-27 0-30 
62 10 181 1100 0-53 0-54 
92 10 202 1110 0-72 0-73 
180 10 280 1100 10 0-97 
269 10 370 1100 10 1-00 
32 3-0 172 1110 0-29 0-30 
62 3-0 184 1110 0-52 0-54 
or 3-0 207 1110 0-70 0-73 
122 340 228 1110 0-82 0-86 
180 3-0 290 1110 10 0-97 
240 3-0 350 1110 10 1-00 
600 2-0 700 1110 1-0 
32 4:5 72 1120 0-30 
g2 4°5 204 1110 0-67 0-73 
151 5 258 1110 0-89 0-93 
210 5 820 1110 1-0 0-99 
300 45 354 1100 10 1-00 
600 45 700 1110 1-0 1-00 
1 100 171 1120 
2-1 100 174 1120 
3-6 100 184 1110 
6-4 100 192 1110 
69 100 195 1120 
10-0 100 205 1120 


*calculated for n = 650 
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(particle diameter ~ 75 — ground 
Sterchamol | brick brick), 
impregnated with a lubricating oil fraction in a 
weight ratio (oil to carrier) of 0-25. Elution 


curves were measured for different concentrations 


300 


(furnace insulation 


of n butane as the feed gas and at varying total 
amounts of feed gas; nitrogen was used as the 
carrier gas. The experiments were performed 
at 0°C and 760 mm exit pressure; the pressure 
drop over the column at the gas rates applied 
was 52mmHg. The data obtained are given 
in Table 1. 

As predicted it is found that the final retention 
volume S,, is fairly constant. Extrapolation 
of the A values to feed 
AS, = ~ 165 ml. Accordingly, nv 
Using eq. (24) n is calculated as 650, With this 


zero volume gives 


~ 1055 ml. 


value the width and maximum height of the 
elution curve can be calculated for the case where 
large feed volumes are used. Comparison of 
calculated and experimental maximum heights 
in Table 1 satisfactory 
Fig. 6 shows the comparison for the width of the 
elution curve. The feed volumes in this figure 
have been plotted as the total amounts of solute 
For the cases of 1-0, 3-0 and 


shows a agreement. 


feed, i.e. n butane. 
WIDTH OF BAND, mi 


l j 
oO 2 6 8 /0 
FEED VOLUME ,mi n-BUTANE 
Fic. 6. The effect of feed volume on band broadening 
in gas — liquid partition chromatography. The curves are 
calculated for n = 650. 


45% of n butane the agreement between theory 
and experiment is excellent*. For pure n butane 
feed the agreement is not so good; strictly 
speaking the theory cannot be applied in this 
case, however, since the volume of the moving 
and immobile phases will not have been constant. 

Fig. 6 that for sharp 
separations (narrow elution bands), the solutes 
the 
concentrated a form as possible. 


shows obtaining 


should be introduced into column in as 
This is easily 
understood, since only in that case is a negligibly 
small fraction of the column length used for 
absorbing the solute, the major part of the column 
being available for the resolution process, 

CONTINUOUS 


THEORY OF THE 


COLUMN 


THE 


General Theory 


The mathematics for the continuous column, 


assuming a linear adsorption or absorption 
equilibrium isotherm and uniform fluid velocity 
over each cross-section, has been given by Lapipus 
and AmuNpson [12]. In this theory longitudinal 
diffusion has also been taken into account. It 
may be remarked here that this longitudinal 
diffusion may partly be of an eddy character, 
i.e. irregularities in flow caused by the packing, 
which are small compared with the size of the 
band and which are statistically distributed 
over the column, may be incorporated in the 
longitudinal diffusion coeflicient. 


As before, we shall consider the case of partition 


chromatography with a linear partition isotherm, 
The 


area (see also Fig. 7) are : 


material balances per unit cross-sectional 


(30) 


The solution of these equations, corresponding 
to the introduction of a concentration pulse 


*The Fig. 6 would coincide if the actual 
volume of diluted feed gas instead of the n butane had 


curves in 


been plotted along the abcissa. 
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VOLUME FRACTION Fy Fy 
CONCENTRATION cy 


O(c, f 
MASS 
WRANSFER| 
d¢ 
LONGITUDINAL CONVECTION 
DIFFUSION 


Fic. 7. Material balance in continuous column. 


Cyt, of sufficiently short duration tf, may be 
derived from the general solution of Lapipus 
and AMUNDSON, 


It reads : 


c ct ut)* at 
ut’)* 
F (t') dt 31 
(31) 
0 
where 
F, F,, (t —t’) 
; 
I, (2, Kt (t t’)) (32) 
| Fy J 


Formulae (31) and (32) may be simplified if 
the column contains a large number of mixing 


stages and transfer units. The number of mixing 


stages for longitudinal diffusion may be defined 


as that number of ideal continuous mixers in 


series which yields the same response to a distur- 


bance in concentration as the given length of a 
column with longitudinal diffusion. From the 
comparison by KLINKENBERG and Ssenrrzer [6], 
it is seen that the number of mixing stages in a 


column of height / is = The height of a mixing 


The 


of transfer, 


stage equals 2D u. number 
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units, or rather the height of a transfer unit, 
is used to characterize the effectiveness of mass 
transfer between two phases. For our case the 
latter may be defined as 
KLINKENBERG and SJENITZER). 
Now it is shown in the Appendix that for 
locations z much larger than both 2D/u and 
F,u/x the equations (31) and (32) reduce to a 
Gaussian function : 


F,u/x (compare 


with 
F ‘ Dz , F? 2 
ei 4 2 2 il 34 
B ' * B K*u (34) 


Since, as already said, D may depend upon 
the flow characteristics of the bed, it is not easy 
to decide to what extent the band broadening 
is caused by longitudinal diffusion on the one 
hand and by mass transfer on the other hand. 
This point will be considered in more detail later. 
The comparison of the Gaussian distribution 
curves as derived from the plate theory and the 
rate theory 


useful relations will be 
derived by comparing the Gaussian distribution 
of the plate theory (17) with the Gaussian 
distribution of the rate theory (33). 

Considering that per unit cross-sectional area 


First a number of 


the effective plate volume is 


v= HF, + au = (35) 
and that 
S = F,ut (36) 
z = nH, (37) 
we find from (17) and (33) that 
uF, 
u 
1+K—! 
h 
— L + 2 (38) 
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where (39) 


is the height of a mixing stage and 


uF, 


(40) 


is the height of a transfer unit. 

Equation (38) shows how the corresponding 
heights for longitudinal mixing (diffusion) and 
mass transfer have to be added to yield the 
H.E.T.P. 

For the the 
Gaussian function combining diffusion and mass 


straightforward derivation of 
transfer mechanisms we are indebted to VAN DER 
(17). This 
from the statistical point of view. 


WAERDEN result is remarkable 
It is evident 
that longitudinal diffusion will lead to a Gaussian 
the 


known that close to the equilibrium (correspond- 


distribution in concentration. It is also 
ing to a large number of transfer units) the fre- 
quently repeated process of transfer of solute 
to the immobile phase and back again may be 
explained in terms of a random walk with result- 
ing Gaussian distribution. 

(33) effects, 


occurring simultaneously, are apparently indepen- 


According to these two when 
dent for statistical purposes, So that the resulting 
Gaussian distribution is found by simply adding 
the 


SJENITZER (6) have computed these component 


compone nt variances. KLINKEN BERG and 
variances and arrived at the same results. 

It is believed that 
applicable in the large majority of practical 
that the 


Amunpbson’s equation to a 


the above conditions are 


cases, SO reduction of LAPIDUS and 


more amenable 


form entails little loss in the field of application. 


Application to the ion exclusion experiments of 


Simpson and Wheaton 


In ion exclusion the ion exchange particles in 
principle act as absorbers for nonionic materials, 
while the ionic material cannot enter, so that the 
ionie material travels at a faster rate through 
the column than the nonionic compounds. 

Simpson and WHEATON [14] carried out a 
series of experiments to determine the H.E.T.P. 


as a function of particle size and flow rate for 


ethylene glycol in water with Dowex 50-X8 
The 


conclusion was that the H.E.T.P. is proportional 


(a cation exchange resin) as absorber. 
to the particle diameter and to the square root 
of the flow rate. The formula, which has been 
used to calculate the H.E.T.P., however, does 
not agree with the plate theory developed in 
the present paper. Since no derivation is given, 
it is impossible to determine how the discrepancy 
arises. 

The main results of Simpson and Wueaton’s 
experiments are summarized in Table 2. 

Table 2. 


Data on ion exclusion of Simpson and WHEATON 


Feed Particle Elution curve 
volume size 


(ml) 


Superficial 
velocity 
Width at 
base 


(ml) 


Maximum 
height 
(frac tion Of Co) 


(mesh) (mm_ sec) 


42 
0-091 32 oF 
0-091 72 
0-091 17 
0-091 20 
0-091 62 
0-091 5 
: 20 
0-091 
0-046 
0-182 
0-365 
O73 


50-100 
50-100 
30-100 
50-100 
50-100 
200-400 
O-70 
30-40 
16-20 
50-100 
50-100 
50-100 
50-100 


2 
5 
5 
5 
5 
5 
5 
5 


The effective column volume nv is 70 ml. 


Fig. 8 shows the influence of the feed volume on 
height and width of the elution curve. It appears 
that 240 and 


reasonable agreement is obtained between experi- 


~~ 
277 a 


with ni 70 ml and n 


ment plate theory. 

For the other experiments, also, the number 
of theoretical plates can be easily determined 
with the aid of Figs. 3 and 4. In this way the 
H.E.T.P. can be determined for every condition 
in two The results are 


independent ways. 


given in Table 3. 
In order to arrive at a quantitative interpre- 
tation of these results we will have to consider 


the mass transfer phenomena more closely, 
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i 


J 
30 40 
FEED VOLUME A, mi 


PARTICLE SIZE 50-/00 MESH 
SUPERFICIAL VELOC/ITY 009/ mm/sec 


J 
/0 20 3O 40 
FEED VOLUME A, mi 


Fic. 8. 
of feed 


Height and width of elution curve as a function 
with the 
results of Simpson and WHEATON. 


volume. Comparison experimental 


Table 3. H.E.T.P. values derived from Simpson 


and Wueraton’s data. 


Number H.E.T.P. 
of plates H 


(mm / sec) n 


Particle 


size 


Superficial 
velocity 
(mesh) (mm) 
200-400 
50-100 
50-70 
30-40 
16-20 


0-091 1000 0-56—- 0-93 18 
0-091 277 20 2-3 10-6 
0-091 : 305 18 - 20 B3-— 8-0 
0-091 155 3-6 a | 10-2 
0-091 3: 49 j 16 


100 0-046 385 5°5-6-6 

100 0-091 : 10-6 
100 0-182 : 14 
100 0-365 2: 3-6 20 
100 0-73 it 33 


The overall mass transfer resistance - may be 
sum of resistances in the 


written as the the 


mobile and the immobile phase 
(41) 


For the mass transfer coefficient in the mobile 
phase the 


ErcGun (18) will be used, which for laminar flow 


correlation formula proposed by 


reads 
25 D,a,? (42) 

6 F, 
D, is the molecular diffusion coefficient in the 
mobile phase, a, is the specific surface area of the 
particles per unit volume of the bed; k, is the 
mass transfer coeflicient per unit surface area 
of the particles, 

For diffusion into a particle the mass transfer 
coeflicient may be defined as 
oc 


ky, (e (43) 


surface 


D,,| 


= Caverage ). 
Or! surface 


A simplified version of the general solution 
for diffusion into a spherical particle can be 
that the 
) does not differ too much 


obtained if it is assumed average 


concentration verage 


from the concentration at the particle surface 
The result then reads : 


(Courtace ). 


a.k 
441 


(44) 
3 d, 


and (44) the overall resistance 


With 


becomes 


(42) 


8 


(45 
22? 


2 

+44) D,a, 

For the case of uniform spherical particles 


(46) 
d 
the 
sum of molecular eddy 
diffusivity FE. In the 
irregular pattern along the particles, D, has to 


The longitudinal diffusivity D will be 
diffusivity D, 


order to 


and 
account for 


be multiplied by a labyrinth factor y which is 
of the order of magnitude of 0-5-1-0. A reasonable 


251 


Co 
/ 
/ 
06 
/ 
WIOTH 
4s 
{mi} 
50 
° 
4 
n= 277 
° 
OL. 50 
956 
20) 7 
4 
/ 
= 
H 
50 
50 
50 
> 


200-400 
3 50-/00 
} 50- 70 MESH 
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25r 


VARYING FLOW RATE 


005 0/0 020 
_mm* sec 


Fic. 9. Determination of longitudinal mixing and particle 
diffusion. 


approach to the eddy diffusivity E is to assume 
E = Aud, (47) 
where A is a packing characterization factor. 
In this way 
D = yD, + Aud, (48) 
Introduction of (45), (46) and (48) into (38) 
yields 


2y (49) 


ud, D, 
where C is a dimensionless factor depending on 
F,, Fy, K and D,/D,, according to: 
F? F,KD, 
75(1—F,P* (1— Dy 
Fy, 


According to (49) the H.E.T.P. expressed in 
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particle diameters is made up of three terms, 
a longitudinal molecular diffusion contribution 
inversely proportional to u times d,, an eddy 
diffusion contribution presumably independent 
of u and d,, and a mass transfer resistance 
contribution proportional to u times d,*. This 
formula suggests that a plot of H/d, versus 


ud, or might decide upon the relative 


ud, 
importance of the three terms. Now by a rough 
estimate it is easily found that longitudinal 
molecular diffusion was of no importance in 
Simpson and Wueraton’s experiments. In Fig. 9, 
therefore, Hd, has been plotted versus F, ud, 
using the values of Table 3. It appears that 
for the experiments with varying flow rate 
(50-100 mesh particles), the linear relationship 
is fairly well confirmed. The experiments with 
varying particle size on the other hand seem to 
follow a completely different curve. The latter 
might, however, be explained by assuming a 
variation of A with the particle diameter as 
indicated by the dashed lines, 

It seems reasonable to expect that it is easier 
to effect a regular packing for large particles than 
for small particles. The values found for A are 
8 for 200-400 mesh, 3 for 50-100 mesh and 
practically zero for 20-40 mesh. 

The particle diffusion coefficient D,,; may be 
calculated from the slope of the straight line 
in Fig. 9. With F, F,, =~ 04 and K 1:3 
it is found that D,, 13 =« 10° m* sec., a 
value which is about 10 times smaller than the 
diffusion coeflicient in the free liquid. 

* In effect, eq. (49) is an example of addition of resistance 
terms. This manner of expression of overall mass transfer 
rates is quite common and a particular analogous case is 
found in the mass transfer in a wetted wall column 
under laminar flow condition. For distillation, WesTHavEerR 
[21] obtained for the H.E.T.P. : 

H Dy, ud 

d wud 192 Dy 
in which d is the diameter of the column. The first term 
in this equation has the same physical significance as 
the first term in eq. (49), whilst the second term is 
essentially the same as the third term in eq. (49). Because 
of the laminar flow conditions and the absence of packing, 
the eddy diffusion parameter A equals zero in WesTHAVER’s 


equation. 
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Application to gas-liquid partition chromatography 

In gas-liquid partition chromatography [19, 20] 
a suitable liquid is retained by the pores of solid 
particles. This liquid absorbs the various com- 
ponents in the flowing gas according to their 
solubility (Henry coefficient), so that a separation 
is obtainable which may be compared with 
extractive distillation. 

In order to describe the diffusion in the liquid 
phase it will be assumed that the liquid within 
the particles is present in the form of a film of 
effective thickness d;. Whether the liquid is 
actually a continuous film or just fills up the 
smaller pores and leaves the wider pores unwetted 
is immaterial for the mathematical description, 
If the latter is the correct concept it will be found 
that the effective film thickness is larger than the 
diameter of the smaller pores. In the first case 
the film thickness should be very small. 

The distribution of fractional volumes of gas, 
liquid, pore space, and particle skeleton in the 
column is schematically indicated in Fig. 10. 


INTERGRANULAR 
SPACE 
GAS -FILLED 
PORES 
| PORE 
| VOLUME 
PARTICLES Liquio 
SKELETON 


Fractional volumes of solid, liquid and gas in 
gas 


Fic. 10. 
liquid partition chromatography. 


The mass transfer coeflicient in the liquid 
phase may again be found from the simplified 
solution of the general equation for diffusion into 


a layer. If d, is the film thickness and a, the 
specific surface of the film, 
1 ,D,,a 


4 d, 


Considering that a,d,;= F,, it follows that 


tes Dy 
2 
4 d; 


(51) 


For the moment it will be assumed that the 
mass transfer resistance in the gas phase may be 
neglected. This assumption will be checked 


later on. Then the general equation (38) becomes 
Fe... KF, ud? 
u 7 A 


The longitudinal diffusivity D will be written 
as the sum of molecular diffusivity yD, and 
eddy diffusivity FE. For the latter approach 
(47) will again be introduced. In this way 


(52) becomes 


) 8 °F, ud? 
u 
1 . 
| il il 


It is seen that, particularly in this case, a 
variation in gas velocity should demonstrate the 
contribution to the H.E.T.P. of gas longitudinal 
eddy and molecular diffusion and liquid molecular 


diffusion. 

We will next consider some experimental 
results obtained with two chromatographic 
columns, both of 180m length and 6mm 


In the first Celite (a diato- 
maceous earth) with an average particle diameter 


internal diameter. 


of about 30, was used as the solid carrier, 
in the second a narrow fraction of ground 
Sterchamol (50-80 mesh; ~ 240). In both 
cases the liquid phase was the same, viz. a lubri- 
cating oil with an average molecular weight of 
312 and a viscosity of 30 cS at 20°C. The weight 
The feed 
sample, consisting of a mixture of normal butane 
(90%) and isobutance (10%), was eluted with 
N, at 30°C and at a pressure close to atmospheric. 
The H.E.T.P. was determined by varying the 


ratio of oil to solid carrier was 0-3. 


sample size and extrapolating to zero sample 
size. The results of such experiments at different 
gas velocities are summarized in Table 4 and 
graphically represented in Figs. 11, 12, 13 
and 14 as a function of gas velocity and reciprocal 
gas velocity. 
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j 
20 3O 


SUPERFICIAL GAS VELOC/TY, mm/sec 


H.E.T.P. versus gas velocity for the celite column. 
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Fic. 12. H.E.T.P. versus reciprocal gas velocity for the celite column. 
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NORMAL BUTANE 


/0 20 30 


SUPERFICIAL GAS VELOC/TY, mm/sec 


Fic. 13. H.E.T.P. versus gas velocity for the sterchamol column. 
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Fic. 14. H.E.T.P. versus reciprocal gas velocity for the sterchamol column. 
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Table 4. Experimental data on H.E.T.P. in 
gas-liquid chromatographic colums. 


Superficial | H.E.T.P. for | H.E.T.P. for 
velocity (normal butane| iso-butane | Solid carrier 
(mm /sec) (mm) 


17 6-0 

20 4-6 

3-4 2-7 

3-5 2-7 

48 2-3 304 

2:3 Celite 

15 18 > 40 
15 20 40-20 52%, 
15 1-7 20-10 » 184% 
15 1-9 10—- 6444% 
23 19 
23 1-9 


1:8 +1 
3-5 
2-5 
2°5 
16 
1-6 
17 240 


1-1 Sterchamol 


(180-300 


For the Celite column the relation between 
H.E.T.P. and gas velocity may be described by 
+ 


9 
H.E.T.P. + O83 + 0-05 u 
u 


for normal butane, and by 


8 
H.E.T.P. = - + 0-8 + 0-08 u 
u 


for isobutane (H.E.T.P. in mm, u in mm_ sec). 


Py 0-019 and 

0-028 for the two components respectively, the 
ratio of the slopes of the straight parts of the 


As in these experiments K 


curves in Fig. 11 according to (53) must be about 
1-5. This appears to be in good agreement with 
the above formulae, where the value of this 
ratio is 1-6. 

The diffusivity D,, in the liquid phase may 


be estimated from its viscosity. This gives 
D,, = 0-3 10° m?/sec. Consequently from 
(53), the known values of the slopes under dis- 
cussion and the A-values mentioned above the 
effective film thickness can be calculated. It is 
found from both curves that d; = 10 p. 

Comparison of the formulae (54) with (53) also 

reveals that 
2yD, ~ 8 10° m?/see. 

This value is of the right order of magnitude as 

y lies between 0-5 and 1, and D, =9 x 10° 

m* ‘see as estimated from the well-known formula 

of GILLILAND [22]. 

The second term in the right hand member of 
(53) is the contribution of eddy diffusion to the 
H.E.T.P. 
the same value for both components. The for- 
mulae (54) do not agree with theory in this 
respect. It must be mentioned, however, that 


It therefore should obviously have 


the values of the constants in (54) are highly 
uncertain, owing to the relatively small number 
of experimental data that were available. The 
value 
2Ad,, 0-5 mm 

seems to be reasonable, from which it follows that 
A = 8. This is in fair agreement with the state- 
ment concerning A on page 282. 

The results of the measurements with normal 
butane in the Sterchamol column may be repre- 
sented by the relation 

H.E.T.P. = 0-05 u + —. (55) 

u 

It will be obvious that the absence of a constant 
in the right hand member is not a definite proof 
that eddy diffusion does not in principle occur. 
Again, with the small number of experimental 
data the constant could not be determined 
accurately. It should not differ much from zero, 
however, so that the eddy diffusion is certainly 
smaller than in the Celite column.* Consequently 


*At a gas velocity of about 17 mm/sec the H.E.T.P. 
of the Sterchamol column just described is about 1-2 mm ; 
the experiments mentioned in the first section of this 
paper show that at identical gas velocity an H.E.T.P. 
value of 2.8 mm is obtained for Sterchamol. This difference 
may well be caused by the difference in homogeneity of 
the packing used in the two cases (50-80 and 50-200 mesh 
respectively). 
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the Sterchamol column is superior to the Celite 
column from an analytical point of view, as 
with less eddy diffusion sharper separations will 
become possible. 

The good agreements between the constant in 
(55) and corresponding constants in the first 
formula of (54) is considered somewhat inci- 
dental. As in the experiments with the Sterchamol 


column K 0-026 the effective film thickness 


Fi, 
in this case is found to be dy 9. With respect 
to the average particle diameter of 240 yu, this 
value seems very reasonable. 

We will finally discuss the question where the 
mass transfer resistance is located. For this 
purpose we will compare the time spent by the 
molecules for travelling to the gas-liquid interface 
with the time needed for diffusing into the liquid 


layer. In the gas phase the molecules have to 


travel over a distance which is roughly one half 


of the pore width. This pore width is of the 
order of magnitude of one fifth of the particle 


diameter, so the path in the gas phase comes into 


the equation as d,. The path in the liquid 
0 


phase is d,, Consequently the ratio of the two 


diffusion times becomes 
d.\? D 
100 \d, dD, 
The resistance in the gas phase may be neg 
lected if this ratio is small in comparison with 
Di, 
100 D, 
10° the two resistances become of equal order of 


unity. As in the present problem 3-3 


magnitude in the case of a film thickness smaller 
than about 0-0001 of the particle diameter. This 
does not the two columns under dis- 
and 


10° for the Celite and Sterchamol columns 


occur in 
cussion. The ratio (56) amounts to 3 
respectively. Apparently the resistance to mass 


transfer is located entirely in the liquid phase. 
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APPENDIX 

The expression (31) will be simplified to a Gaussian 
distribution function by introducing a few approximations, 
which are essentially due to Van per Warrpen [17]. 

In the first place the first term of the asymptotic 
expansion of the Bessel function J, will be used. This 
involves that the argument of J, is sufficiently large, a 
condition which will be verified afterwards. Then (32) 
reduces with 
1 Fy 
K Fy 


aK?F, t—t/B 
Fy V 


1/8 


Kt’ 
24 Fy, Fy, (t 


2) 
ve | 


Now in the integral (31) two exponential functions 


exp 


occur. Their maximum values are obtained for (’ 


and f Bt respectively. We will next assume that the 
effects of longitudinal diffusion and of the deviation from 
equilibrium are small. This 


means that all retention 


times ¢ are close to the mean { 2° which is that prevailing 
u 


in ideal linear chromatography. Consequently the two 


maxima of the individual exponential functions ¢ 


and ¢ At are close together, so that the main contribu- 
tion to the integral depending on the product of those 
exponential functions will also be located in this small 


region where ( = ~ = Bt. It follows that it is permissible 
u 


to replace t/ by — or Bt, except in the numerators of the 
u 


exponents. In this way F (t’) may be further simplified to 


1 Fy, 
exp 
| 


It also follows that the first term of (31) may be com- 
pletely neglected, as long as z ~ 8 ut is markedly different 
from ut, i.e. if there is a sufficient retention of the solute 
in the column. 


xK?F, t’\ 2) 
Fit) (: 
B 


Now the expression (31) may be simplified to 


(z/u-—8ty* 
2(o1* + os") 


where 


to 
956 
- 
4 
» 
= 
Bly e e dt 
0 
2 2 e 
V 2m + 03°) 
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K*u 


2 2 9 a2 
a) 2 2, 


The 3 assumptions used in this derivation may now be 


specified as 


a* Kt (t —t') x2 
2 2 F u? 
P,P yy 


«\2 ‘ 
2D 
< or <1 
u uz 


2u 


The first condition requires that the height of a transfer 


z 2 
< = or 
u 


F yu 
unit is much smaller than the height z of the bed. 
x 
According to the second condition the height of a longi- 


tudinal mixing stage — 


the height of the bed. The last condition is always fulfilled 
if the first one is valid, because 


should be much smaller than 


BF yy 1 
KF, (1 LK 

It follows that the simplified expression may be used 
if the column contains a large number of mixing stages 
and transfer units, or, in other words, a large number of 
theoretical plates. 
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Abstract—The heat flux from an electrically heated horizontal platinum wire to water, to 
methylethylketone and to an aqueous mixture containing 41°, wt. of methylethylketone was 
determined throughout the ranges of convection and nucleate boiling, both at atmospheric pressure 
and at a pressure of 42 em Hg. 

A general shift of the curves to lower heat transfer occurred with decreasing pressure, as a 
consequence of the increasing average size of vapour bubbles, generated on the heating surface. 
The maximum in nucleate boiling heat flux to the 41%, mixture exceeded considerably the 
corresponding value in water at both pressures investigated. 

Maximum heat flux in nucleate boiling has been determined systematically as a function of 
concentration in a range of concentrations with water in excess, for mixtures of water and 
methylethylketone, acetone and the lower aliphatic alcohols respectively at several subatmo- 
spheric pressures. 

In all mixtures a maximum value of the maximum heat flux occurred at a certain low 
concentration of organic component, which was approximately independent on pressure. The 
ratio of these maxima to the maximum value in water af the same pressure decreased with 
decreasing pressure. 


Résumé— Pour déterminer le coeflicient de transmission thermique d'eau, de méthyl-éthyl- 
cétone et dun mélange contenant 41 


auteur a utilisé un fil horizontal de platine servant d’élément chauffant. Il a mésuré les diverses 


de m.-é-c, bouillant sous les pressions de 76 et 42 em Hg, 


valeurs de transfert jusqu’au maximum “ nucleate boiling.” 
Pour une réduction de la pression le transfert décroit 4 cause d'une augmentation de la 
* nucleate boiling 


dimension des bulles de vapeur générées au fil chauffant. Le maximum 


du mélange & 41°, dépasse considérablement la valeur correspondante pour eau sous les deux 
pressions. 

L.auteur a mésureé les valeurs du flux de chaleur au maximum “ nucleate boiling ~ en fonction 
du concentration pour les mélanges d'eau et de m.-é.-c., d’acétone et des premiers alcools aliphati- 
ques, sous diverses pressions subatmosphériques. L'auteur resteint a étude des mélanges 
contenant un excés d'eau. 

Pour tous ces mélanges un maximum élevé du flux de chaleur a été trouvé pour une petite 
concentration du composant organique, qui est approximativement indépendant de la pression, 
La proportion de ces maxima au maximum correspondant d'eau décroit pour une réduction de 


le pression. 


Si 
The rate of heat flow q from a heated wire to a 
T. the 


difference in temperature between the heating 


RVEY 
boiling liquid is a function of 0 t 


surface and the bulk liquid, The curves represent- 


a(@) 
flux 


ing the heat where A denotes the 


of three 


convection, 


surface, 


of the 
different parts : 


heating consist 
the of 
which no vapour bubbles are observed, the region 


area 
region in 
of nucleate boiling, in which vapour bubbles are 


generated on special spots of the surface, and the 


of film boiling, in which a gradually 


increasing area of the surface is covered witha 


region 


coherent layer of In the regions of 


convection and nucleate boiling the values of 


vapour, 


@ are small or moderate, in contradistinction to 
the region of film boiling, where higher values 
of @ occur, up to the melting point of the wire 
(3, 6). 


Experiments, had been carried out 


which 


previously at atmospheric pressure [5], showed 


that the maximum of nucleate boiling in a mixture, 


containing 42°, wt. of methylethylketone in 


200 


q 
f 
x 
= 
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water, exceeded considerably the maxima in the 


pure components. This experiment was now 
repeated more accurately*, while investigations 
at a lower pressure were made at the same time. 
heat 


atmospheric pressure had also been determined 


Maximum flux in nucleate boiling at 


previously as a function of concentration for 


mixtures of water with acetone, 
methylethylketone the 
respectively [6]. In all 


value of the maximum heat flux occurred at a 


binary 
lower alcohols 
systems a maximum 
certain low concentration of the organic com- 


ponent. Similar measurements at several 


subatmospheric pressures are carried out at 
present for these mixtures in a range of con- 


centrations with water in excess. 


METHOD AND APPARATUS 


The method described by McApams [3] has 
in principle been used to determine the heat 
transfer curves up to the maximum of nucleate 
boiling for water, for methylethylketone (MEK) 
and for a mixture containing 4-1°,, wt. of MEK 
in water, both at atmospheric and at lower 
pressures, 


For that 
passed through a horizontal platinum wire, with a dia- 


purpose a gradually increasing d.c. was 
meter of 0-02 em, of which a central portion was used 
as a test section and provided with thin platinum potential 
taps. This section served as the heating surface, and at the 
same time as a resistance thermometer [3, 5]. The test 
sections had a length of approximately 5 em, and the 
diameter of the potential taps was 0-005 em. The wires 
were annealed at red heat (1,000°C) for 15 min. 


The 


maximum nucleate boiling heat flux was investi- 


influence of concentration upon the 
gated systematically at several subatmospheric 
pressures for aqueous binary mixtures of MEK, 
acetone and lower aliphatic alcohols. 
Restriction was 


water in excess, and in this limited 


some 
made to mixtures containing 
range of 
concentrations results are compared with data 


obtained previously [6] at atmospheric pressure. 


*Owing to a calibration error, the temperatures given 
in the previous papers [5-Figs. 2, 3 and 4) and [6-Figs. 
2 and 3] are too high by approximately 10°C. In conse- 
quence the values for the coefficient of heat transfer 
h [5-Fig. 4] and h [6-Fig. 3] respectively, are also 
erroneous. 


max 
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In these experiments a.c. was used for heating, as the 
temperature difference between wire and bulk liquid 
was not measured. The wires burned out in some cases, 


when relatively high maxima were obtained. 

The boiling vessel was in principle similar 
to that described previously [5]. It consisted of 
The 


nickel-coated brass cover ring and bottom plate 


a glass cylinder with ground end faces, 


were fastened to the cylinder by means of draw 
bars, and neoprene rings were used for packing. 
A pertinax cover was provided with four nickel- 
The 
while the potential leads 
attached to the 
were placed at a suitable distance from the wire. 
A total 


cover plate, and the boiling vessel was only 


coated brass bars. wire was stretched 


between two bars, 


were remaining bars, which 


reflux condenser was mounted on the 
partly filled with liquid; it was heated from 
below by a bunsen flame to the boiling tempera- 
ture desired. The combined action of a water 
jet pump and a Cartesian manostat ensured a 
constant vacuum. Uniform boiling was main- 
tained, when the experiments were carried out, 
as heat losses were compensated by the flame. 
In all cases the boiling temperature could be 
kept constant within 0-5 °C. 


The MEK and the 


1-pentanol were chemically pure, the latter was 


Distilled water was used. 


grade, the 1-butanol had a 
115-117 C, and the 


organic liquids were of commercial grade. Boiling 


the furfural free 


boiling range from other 
and equilibrium data, pertaining to the system 
water — MEK at atmospheric and subatmospheric 


pressures, have been taken from literature [2, 4]. 


, ( 
Water. The heat flux * to boiling water was 


as a function of the temperature 
Results at 
The 


figures at the curve denote the number of nuclei 


determined 
difference 90 to serve as a reference. 


atmospheric pressure are shown in Fig. 1. 


at which vapour bubbles were generated on the 
The test 
The 


increased gradually in the region of nucleate 


heating surface. area was 0-31 cm? 


approximately. number of active nuclei 


boiling, until film boiling started. The average 


diameter of the bubbles leaving the wire amounted 


: 
OL. 
956 
4 
4 
ay. 
RESULTS 
i 


S. J. D. VAN STRALEN 


to 02cm. The reproducibility of the measure- 
ments at atmospheric pressure was satisfactory, 
as follows from the data obtained with two 
different heating wires. 

The heat flux to water at pressures of 42, 20 


76 erm Hg 


—e fy 


Fic. 1. Heat flux for convection and nucleate boiling to 
water as a function of temperature difference between 
heating surface and liquid at atmospheric pressure and at 
pressures of 42, 20 and 10cm Hg. The figures at the 
eurves denote the number of nuclei generating vapour 
bubbles on a platinum heating wire with a length of 
approximately 5cem and a diameter of 0-02 cm. 
Wire 1 : 76cm Hg (V), 10 em Hg ( @) ; wire 2 : 20 em Hg 
42cm Hg (OC); wire 3: Técm Hg (YF). 


of 
max. (cal sec-' 


= 


20 


and 10cm Hg had also been determined, in the 
latter case with one of the wires, which was used 
already at atmospheric pressure. In this case 
only two active nuclei were present, generating 
vapour bubbles with relatively large diameters 
of 1-lem approximately (Fig. 1). The slope 
in the region of nucleate boiling tallied fairly 
well with that in the natural convection part 
of this curve. 

The larger size of vapour bubbles was respon- 
sible for the lower maximum of nucleate boiling 
at a lower pressure. This is clearly illustrated 
by the shape of the curve, showing maximum 
nucleate boiling to water as a function of pressure 
p (Fig. 2). The maximum heat flux increased 
continuously with increasing pressure, while 
the average size of the bubbles decreased at the 
same time. The maximum of nucleate boiling 
approached to the relatively small value of 8 cal 
at pressures below 10cm Hg. This 
phenomenon was caused by the fact, that film 
boiling occurred immediately at first appearance 
of only one vapour bubble, with relatively large 
diameter, on the wire in this case. 

Methylethylketone and acetone. In Fig. 3 the 
heat flux to water, and to 4-1 and 100°, wt. of 
MEK is shown as a function of the temperature 
difference between heating surface and bulk 


40 


—= pi(cm Hg) 


Fic. 2. Maximum of nucleate boiling to water as a function of pressure. 
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Heat flux for convection and nucleate boiling 
and 100° wt. of methylethylketone 


Fic. 3. 
to water and to 41°, 
as a function of temperature difference between heating 
surface and liquid at atmospheric pressure and at a pressure 
of 42 em Hg. The figures at the curves denote the number 
of nuclei generating vapour bubbles on a platinum heating 


wire with a length of approximately 5 cm and a diameter 


of O-O2 em. Water(...... ) wires 1 and 3: 76cm Hg; 
wire 2: 42cemHg. 41% wt. of MEK ({ ) 
wire 4: Hg (7); wire 5: 42cm Hg (YW). MEK 
( ) wire 4: 76 cm Hg (QO) ; wire 6 : 42 em lig (@). 


liquid, at atmospheric pressure and at a pressure 
The maximum heat 
t-1°, wt. of MEK 


the corresponding value for water 


of 42cm Hg respectively. 
flux in nucleate boiling to 


exceeded 


Heat transfer to boiling binary liquid mixtures at atmospheric and subatmospheric pressures 


Fic. 4. 
as a function of composition at 76 cm Hg ( 


Wom Hg and 29cm Hg (@ ). 
carried out with the same wire are represented by the 
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both 
This behaviour could be explained by the fact, 


considerably at pressures investigated. 
that the average size of vapour bubbles, leaving 
the heating surface, was smaller in this mixture 
than in water [6]. 

A general shift of the curves to lower heat 
transfer at the same temperature difference, 
and lower maximum values in nucleate boiling 
occurred with decreasing pressure. 

The maximum of nucleate boiling was deter- 
mined systematically as a function of concentra- 
water—-MEK in a 
limited concentration range with water in excess, 
10 and 29cm Hg (Fig. 4). 


The data, which were obtained at atmospheric 


tion for the binary system 


at pressures of 76. 


pressure, have been published previously [6] and 


values. All 


show a maximum value at approximately the 


are included as_ reference curves 


same concentration, whereas the absolute value 
of this 


decreasing pressure, 


maximum decreased gradually with 


A similar behaviour was found with the system 


| 
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Measurements 


Water-methylethylketone. 


same figures. 
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water-acetone, which was investigated at pressures 
of 76 and 30cm Hg (Fig. 5). 
maxima of the curves also occurred at the same 


In this case the 


concentration. 
These results are in fair agreement with the 
theory [6], as the ratio of the equilibrium constant 


30cm Hg 


"5 20 
Wi % Acetone 


Fic. 5. Water-acetone. Maximum heat flux as a function 

of composition at 76cm Hg (O © @ (—)) and 30cm Hg 

(B \ A @). Measurements carried out with the same 
wire are represented by the same figures. 


£0 


20 
—e W! % Ethanol 


W ater-ethanol. Maximum heat flux as a function 
29 em Hg 


Fic. 6. 
of composition at 76em Hg (VY Vv ) 


(@O). 


and 
Measurements carried out with the same wire 
are represented by the same figures. 


£20 
—e W! Y% !- Propanol 


Water-1-propanol. Maximum heat flux as a func- 


Measurements 


Fig. 7. 

tion of composition at TéemHg 

Hg (O Y) and 10cm Hg (OW). 

carried out with the same wire are represented by the 
same figures. 
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—e Wt % Butanol 
Maximum heat flux as a function 
(® 4), 32 em Hg 
A) and 10cm Hg (O @). 


wire are re- 


Fic. 8. Water-1-butanol. 
of composition at 76 em Hg 


(a V SO). 19 em Hg ( 


Measurements carried out with the same 


presented by the same figures. 
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Heat transfer to boiling binary liquid mixtures at atmospheric and subatmospheric pressures 


for the most volatile component to the equilibrium 


constant of water, at a fixed concentration, is 
only slightly dependent on the pressure in the 
10-76cem Hg. This 
literature [4] for both systems. 


The lower aliphatic alcohols. The maximum of 


range from follows from 


nucleate boiling to the binary systems water- 
ethanol (Fig. 6), water-1-propanol ( Fig.7), water-1- 
butanol (Fig. 8) and water-1-pentanol (Fig. 9) 
was determined as a function of concentration 
in a limited range of concentrations with water 
These 


subatmosphe ric 


were made at 


Results at 


previously 


In EXCESS. investigations 


several pressures. 


atmospheric pressure, which were 
obtained [G), are included as reference values. 


maximum 


these systems a maximum value in 


nucleate boiling heat flux occurred 


7 8 


Wit % 1-Pentanol 


flux as a 


27 cm Hg 


Maximum heat 
TéemuHg (= @ 
© @). 


with the 


Fic. %. 

function of 

(OF). 


Water-1-pentanol, 
composition at 


Measurements carried out 


and 
same wire 


are represented by the same figures. 


at a certain low concentration of organic com- 
ponent, which was approximately independent 
on pressure. The absolute values of the maxima 
decreased gradually with decreasing pressure. 

At a pressure of 10cm Hg no maximum were 
the curves (Figs. 7 and 8), as the 
size of vapour bubbles had 
the that the heating 


immediately into film boiling if only one bubble 


present in 
been increased to 


extent wires passed 


appe ared. 


maximum heat thax 


Fic. 10. Ratio of 


mirtares to corresponding value for 


cater as a Ttunctr 


pressure 


heat flux to that 


the 


Ratio of 
Ne 


nucleate 


the 


boiling 


absolute values of 
heat 
with pressure, but even their ratio to the maximum 
th This effect 
10 tor ted. 


The curves pe rtaining to both the alcohols and 


only maxima 


flux inereased eradui 


value in water at san 


pre ssure, 


is shown in Fi all svstems Investig: 


the ketones have a similar shape, and the effect 
clearly depends primarily on the number of carbon 


atoms In the molecules concerned. 
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